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.^The  present  investigation,  which  was  primarily  of  an  experimental  nature, 
was  undertaken  to  define  the  flow  field  produced  when  an  underexpanded  jet  Is 
exhausted  from  a conical  nozzle  (rocket  exhaust)  into  a cylindrical  tube  (tube 
launcher).  The  expcrimrntal  data  Included  the  reservoir  conditions  for  the 
nozzle,  the  static  pressure  distribution  along  the  cylinder  wall,  pit'.t  probe 
su-veys  at  the  cylinder  exit,  and  photographs  of  the  flow  exhausting  from  the 
tube.  Of  principal  concern  were  the  effects  on  the  flow  field  of  tube  length, 
of  reservoir  conditions,  end  of  test  time.  Although  room- temperature  nitrogen 
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DEFINITION  OF  SYMBOLS 


Symbol  Definition 

A Area  of  spherical  surface 

D Internal  diameter  of  tube 

K Mach  number 

m Number  of  Increments  or  iterations 

P Pressure 

R Spherical  radius 

r,  z Cylindrical  coordinates 
S Entropy 

T Temperature 

8 Angle  of  the  impingement  shock  with  respect  to  the  tube  axis 

Y Ratio  of  specific  heats 

A9  An  increment  of  the  nozzle  half  angle 

9 Angle  of  the  flow  with  respect  to  the  tube  axis 

9 Conical  nozzle  half  angle 

n 

U Mach  angle 

V Prandtl-Meyer  angle 
Superscript 

* Evaluated  at  the  nozzle  throat 

Subscript 

a Evaluated  at  atmospheric  conditions 

b Evaluated  In  the  dead  air  region 
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DEFINITION  OF  SYMBOLS  (Concluded) 

• 

Evaluated  at  stagnation  conditions 
Evaluated  at  point  p 
Evaluated  at  the  tube  wall 
Evaluated  along  the  tube  wall 
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CHAPTER  I 


INTRODUCTION 

A primary  problea  aasoclatcH  with  the  structural  design  of 
lightweight  cylindrical  tube  rocket  launchers  is  the  prediction  of  the 
maximum  pressures  or  forces  exerted  on  the  tube  wall  by  the  flow 
exhausting  from  the  rocket  nozxle.  An  understanding  of  the  rocket 
exhaust  flow  field  in  the  launcher  is  a prerequisite  to  the  determina- 
tion of  the  pressure  distribution  along  the  tube  wall  and  to  the 
optimum  launcher  structural  design.  The  current  demand  for  more  power- 
ful rockets  coupled  with  the  requirement  for  lightweight  launchers  has 
only  recently  emphasizjd  the  need  for  a knowledge  of  the  flow  field 
inside  the  launcher.  Up  to  now,  the  maximum  pressure  in  the  tube  was 
assumed  to  be  that  at  the  nouzle  exit.*  The  structural  integrity  of 
launcher  was  usually  demonstrated  by  firing  the  rocket  from  the  launcher. 
Some  experimental  information  was  obtained.  However,  it  was  not  used 
to  gain  an  understanding  of  the  exhaust  flow  but  to  ascertain  if  the 
launcher  could  withstand  the  pressures.  Therefore,  the  research 
reported  In  this  paper  has  been  directed  to  the  understanding  of  the 
rocket  exhaust  flew  In  the  launcher  nr.  In  other  words,  of  an  axlsym- 
metrlc,  underexpanded  Jet  exhausting  Into  a cylindrical  tube. 

A literature  survey  which  was  conducted  through  NASA's 
Scientific  and  Technical  Information  Facility  and  the  Defense 
Documentation  Center  failed  to  provide  specific  Information  regarding 

* US  Army  Missile  Coinnand  Report  No.  RL-T31-71-9,  dated  30  July  1971. 

5 


FRECEDItO  PAOS 


,r»T  f iiieD 


th*  rc.;k*t  eif--*-*:  !i*ii  * CTllruirSc*!  fJt*  I*^=:h*r.  Hc%r«-irer, 

tr.fcrajit  Ic^a  prrtaialr.j  to  ctr-#r  w*»  vfhlch  v«s 

to  thik  a^^cific  a;;ltcatlca.  A flew  ac-^el  b-LscNi  u^ca  t^v« 
ftndlr.gs  cf  ta#  llt*rat’-r«  **4r;h  w*»  4evi»<-i  to  th«er«tl c*lly  *p?roxt- 
a-»t#  th«  flcv  vithlR  th«  t-jfce.  A co*;fjt*r  progrAa  vhich  utillxed  the 
aethod  of  characterttt Ics  for  steadT.  aztayanet ric . rctatlcAal  flow  and 
which  was  forxulated  frea  t" at  flrw  r«odcl  wa*  developed  to  obtain  the 
pressure  distribution  along  the  tube  wall  and  to  gain  an  Insight  Into 
the  characteristics  of  the  flow.  For  conditions  behind  the  first  shock, 
an  approxlration  taethod  was  Incorporated  Into  the  progran  to  elisinate 
the  necessity  of  dcterelnlng  the  entropy  variation  in  the  noresal  direc- 
tion along  each  streasllne.  The  base  pressure  obtained  fros  experi- 
ments and  from  an  existing  theoretical  laethod  based  upon  Korst's  flow 
t»del  was  used  as  an  input  into  the  corputer  prograa. 

The  theoretical  investigation  was  performed  to  augment  the 
experimental  investigation  which  constituted  the  aajcr  portion  of  this 
research.  The  experimental  Investigation  consisted  pritaarily  of  tests 
employing  cold  gas  flows,  although  some  tests  employing  hot  gas  flows 
were  performed.  A cold  gas  facility  using  dry  nitrogen  was  desired 
and  fabricated  to  develop  a supersonic  underexpanded  Jet.  Double  base 
solid  propellant*  rocket  motors  produced  similar  jets  for  the  hot  gas 
tests,  ihese  Jets  were  exhausted  Into  instrumented  tubes  of  different 


* A double  base  solid  propellant  is  ermposed  mainly  of  a colloid  of 
nitroglycerin  and  nitrocellulose. 
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lengths  to  simulate  the  various  positions  of  the  rocket  in  the  launcher 
during  launch.  High  frequency  response  pressure  transducers  were  used 
to  sense  static  wall  pressures  along  the  tubes  and  to  obtain  pitot 
pressures  within  the  cold  gas  flows.  High  speed  pictures  were  taken  of 
the  flow  exiting  from  the  tube  to  aid  in  the  analysis  of  the  flow 
within  the  tube. 

Finally,  the  results  from  the  present  experimental 
Investigation  and  from  similar  investigations  available  in  the  litera- 
ture are  correlated  r/lth  the  results  of  the  theoretical  Invr jtigatlon. 
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CHAPTER  II 
LITERATURE  REVIEW 

In  1940,  Nusselt  (1]  conducted  experiments  with  supersonic 
gas  (air)  discharging  into  an  enlarged  circular  cross-section  tube. 
When  there  is  an  abrupt  increase  In  flow  cross-scc tional  area  as  shown 
in  Figure  1,  an  isolated  subsonic  region  (dead  air  region)  will  be 
created  at  the  corner.  Tf  viscous  effects  are  neglected,  the  pressure 
in  this  region,  Pb  (referred  to  as  base  pressure  by  Wicks  [2]  is 
impressed  upon  the  jet  boundary  and  is,  in  fact,  the  jet  boundary 
pressure.  In  reality,  the  conditions  in  the  dead  air  region  are 
determined  primal lly  by  viscous  effects  along  the  jet  boundary, 

Nusselt  concluded  that  Che  base  pressure  could  be  equal  to,  less  than, 
or  greater  than  the  nozzle  exit  pressure,  P^,  depending  on  the  flow 


Figure  1.  Flow  Field  for  an  Axisymmetric , Underexpanded, 
Supersonic  Jet  Exhausting  into  a Cylindrical  Tube 
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total  pressure.  The  description  of  the  flow  fields  for  which  these 


conditions  exist  was  not  given. 

In  1955,  Fabri  and  Slestrunck  (31  published  a comprehensive 
study  of  a supersonic  stream  entering  an  abrupt  circular  enlargement. 
The  type  of  flow  within  the  enlarged  section  or  tube  varied  with. an 
increase  in  total  pressure  from  mixed  flow  (subsonic  and  supersonic  in 
the  tube)  with  'low  separation  in  the  nczzle  to  fully  developed 
supersonic  flow  throughout  the  tube.  As  did  Nusselt,  they  used  the 
base  pressures  and  total  pressures  to  define  Che  flows  .~nd,  by  the 
use  of  shadowgraphs,  were  able  to  associate  these  pressures  with 
visual  obsen'atlons  of  the  flow  fields.  Shadowgraphs  obtained  by 
Fabri  and  Slestrunck  (3]  of  the  flow  field  with  the  nozzle  exhaust 
filling  the  tube  are  reproduced  in  Figure  2.  A coin>lete  sequence 
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of  the  flows  are  shown  in  references  [3],  [4],  and  (5).  The  base 
pressure/total  pressure  relation  corresponding  to  the  shadowgraphs  is 
shovm  in  Figure  3.  Figure  3 shows  that  the  base  pressure  decreases  as 
the  total  pressure  increases  in  the  mixed  flow  regime  and  increases  as 
the  total  pressure  increases  in  the  supersonic  regime. 


The  tra.isir.ion  from  mixed  flow  to  fully  supersonic  flow 
begins  immediately  when  the  Jet  fills  the  tube.  The  shadowgraph 

pictured  in  Figure  2. a shows  the  Jet  shear  layer  in  contact  with  the 
wall.  A series  of  lambda  shocks  form  in  the  tube  reducing  the 
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supersonic  flow  to  subscnic.  The  subsoric  shear  layer  along  the  tube 
wall  provides  a relatively  large  means  of  communication  for  the 
downstream  conditions  back  to  the  dead  air  region.  The  corresponding 
transition  base  pressure  is  indicated  in  Figure  3.  Transition  begins 
at  a base  pressure  which  is  lower  than  the  nozzle  exit  pressure 
indicating  that  the  nozzle  exhaust  flow  is  expanding  from  the  nozzle 
exit  diameter  toward  the  larger  tube  diameter. 

An  increase  in  total  pressure  causes  the  base  pressure  to 
decrease.  The  gases  in  the  jet  boundary  expand  and  impact  the  wall. 

An  oblique  shock  wave  forms  (Figure  2.b),  but  the  shock  pattern  breaks 
down,  indicating  that  supersonic  flow  does  not  exist  throughout  the 
tube.  Once  the  total  pressure  is  sufficiently  increased,  the  flow  in 
the  tube  is  supersonic  and  Intersecting  oblique  shock  patterns  form 
throughout  the  tube  (Figure  2.c).  Note  that  the  upstream  flow 
characteristics  in  the  supersonic  flow  should  now  become  essentially 
independent  of  downstream  conditions.  Crocco  [5]  noted: 

Observe  that,  as  the  pictures  (refer  to  Figure  2} 
show  clearly,  the  shocks  do  not  reach  the  wall  because  of  the 
presence  of  subsonic  flows  in  the  boundary  layer  and  in  the 
zone  where  the  Jet  mixes  with  the  gases  recirculating  in  the 
dead  water  (air)  region.  However,  these  subsonic  regions  are 
thin,  and  if  their  longitudinal  extent  is  sufficient,  they 
cannot  prevent  the  Isolatlun  of  the  dead  water  region  from 
the  downstream  effects. 

The  comparison  of  theory  with  experimental  data  presented  in 
Figure  3 indicates  that  the  one-dimensional  mixed  flow  theory  of  Fabri 
and  Siestrunck  is  adequate  for  the  mixed  flow  regime.  However,  one- 
dimensional theory  no  longer  applies  with  the  advent  of  a complicated 
axisymmetric  supersonic  flow.  Also,  the  onset  of  supersonic  flow  cannot 
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be  predicted  with  one-dimensional  theory.  According  to  Fabri  and 
Siestrunck,  the  base  pressure/total  pressure  relation  for  fully 
developed  supersonic  flow  is  Pjj/Pq  * constant.  The  value  of  this 
constant  must  be  obtained  either  from  experiment  or  from  a theory 
which  involves  the  consideration  of  the  shear  layers  bounding  the 
dead  air  region  along  with  the  characteristics  of  the  supersonic 
flow. 

Fabri  and  Siestrunck  recorded  the  tube  wall  pressures  and 
presented  the  ratio  a function  of  tube  length  for  the  super- 

sonic regime  (Figure  4).  The  results  of  their  axlsymmetric  method  of 
characteristics  study  are  also  shown  in  Figure  4.  The  characteristics 
net,  which  is  not  presented,  was  shown  as  a mirror  image  of  a shadow- 
graph similar  to  Figure  2.c  and  included  the  shock  emanating  from  the 
Jet  boundary-wall  intersection  and  its  reflection  from  the  centerline 
back  to  the  wall.  No  information  was  provided  regarding  the  details  of 
solution  of  the  complex  axlsymmetric,  rotational  flow  with  multiple 
shocks . I 

An  interesting  practical  application  of  the  type  of  flow 

I 

1 

being  investigated  is  found  in  wind  tunnels  with  a free  Jet  test  sec- 
tion. The  particular  application  was  treated  theoretically 
(one-dimensional  theory)  ty  Hermann  [6].  As  can  be  noted  from  later 
publications  of  Fabri  and jsiestrunck  [4,  7],  another  associated 
practical  application  is  the  supersonic  ejector.  In  1958,  Coethert  [81 
proposed  that  an  .ejector  which  utilized  the  test  engine  exhaust  (xocket 
nozzle  exhaust)  could  be  used  to  provide  high  altitude  simulation 
pressures  in  a test  cell.  This  triggered  a flurry  of  activity  by  the 
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U.S.  Air  Force  at  Arnold*Alr  Force  Scation,  Tenncsaee,  and  Edwards  Air 
Force  Base,  California  [9  - 18]. 


Figure  4.  Versus  Tube  Length  for 

Fully  Developed  Supersonic  Flow  [3] 

The  ejector-diffuser  application  is  based  upon  the  relatively 
low  base  pressure  obtained  when  exhausting  a supersonic  stream  into  an 
enlarged  section.  The  minimum  base  pressure  occurs  at  the  beginning  of 
the  supersonic  regime  (Figure  3).  Most  of  the  work  in  this  field  has 
revolved  around  understanding  the  base  pressure  phenomenon  in  order  to 
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obtain  the  minlmutn  base  pressure  possible.  Consequently,  the  character 
izatlon  of  the  supersonic  tube  flow  has  beer  neglected,  especially  at 
higher  total  pressures  than  those  necessary  to  produce  nilnimuai  base 
pressure. 

The  adaptation  of  Korst's  base  pressure  flow  model  [19]  to 
axlsymmetrlc  ejectors  without  Induced  flow  by  Bauer  [16,  17]  has 
evolved  as  the  best  rethod  of  predicting  base  pressure.  This  base 
pressure  prediction  method  is  also  applicable  to  the  specific  problem 
posed  in  this  work.  Other  pertinent  theories,  i.e.,  those  of  Crocco 
and  Lees  [20],  Zumwalt  [21],  Nash  [21],  Greenwood  [22],  and  Lamb  [23] 
may  be  adaptable  to  this  specific  problem.  Bauer  [17]  and  German  et  al 
[18],  proposed  a method  for  the  determination  of  the  tube  wall  pressure 
at  the  jet  boundary-wall  intersection.  This  pressure,  termed  by  Bauer 
as  the  peak  recompresslon  static  pressure,  is  also  based  on  Korst's 
theory.  This  prediction  method  agreed  favorably  with  their  experi- 


mental data. 


CHAPTER  III 


ANALYTICAL  INVESTIGATION 
Flow  Model 

A highly  underexpanded  supersonic,  short  duration  Jet 
exhausting  from  a conical  nozzle  Into  a short  tube  with  an  Inside 
diameter  slightly  larger  than  the  nozzle  exit  characterizes  the  flow- 
field  of  tube  launched,  folding  fin  rockets.  A simple  flow  model  has 
been  devised  to  describe  such  a flow.  The  model  is  based  upon  the 
fully  developed  supersonic  flow  as  defined  by  Fabrl  and  Slestrunck. 

The  operating  rocket  motor  chamber  pressures  vary  according 
to  the  burning  characteristics  of  the  grain.  The  maximum  pressures 
range  from  250  to  350  (P»/P  },  depending  on  the  rocket.  As  can  be 

V A 

noted  from  Figure  3,  these  total  pressures  are  far  In  excess  of  that 
required  to  produce  supersonic  flow.  Foster  [11]  experimentally 
investigated  flows  with  maximum  Fq/P^  " ^0.  The  large  total  pressures 
coupled  with  the  geometry  should  produce  fully  developed  supersonic 
flow. 

In  the  present  model  steady  flow  is  assumed.  However,  in 
reality,  the  flow  is  not  steady  during  rocket  motor  Ignition  or  chamber 
pressure  talloff.  Therefore,  the  analysis  will  concentrate  on  those 
data  obtained  during  the  quaoi-steady  flow  which  exists  during  the 
Interval  when  the  stagnation  pressure  Is  at  or  near  a maximum.  The 


flow  field  "during  this  time  interval  is  of  Interest  from  an  engineering 
standpoint,  because  the  stresses  in  the  tube  wall  are  at  a maximum. 

Heat  transfer  and  real  gas  effects  (including  possible  chemical  reac- 
tions in  the  tube)  are  neglected.  However,  they  should  be  included  if 
a more  rigorous  analysis  of  the  hot  gas  flow  is  desired.  The  flow  is 
assumed  to  be  inviscid,  although  a theory  wtiich  includes  the  viscous 


nature  of  the  gas  at  the  jet  boundary  is  used  to  predict  base  pressures 
Figure  5 is  a schematic  of  the  flow  model. 


viscoos  jrr  inviscid  jet 


Figure  5.  Flow  Model  for  Axisymmetric,  Underexpanded , 
Supersonic  Jet  Exhausting  into  a Cylindrical  Tube 


16 


The  radial  flow  region  bounded  by  the  nozzle  exit  spherical 
radius,  R^,  Prandtl-Meyer  expansion  fan,  and  Impingement  shock  is 
termed  the  conical  flow  region.  The  characteristics  of  the  flow  In 
this  region  are  determined  from  geometry  and  the  stream  tube  relation, 

R 

-or  " f(Y.M  ) . (1) 

The  axially  symmetric,  Isentroplc  flow  region  is  bounded  by 
the  Prandtl-Heyer  expansion  fan,  Inviscld  Jet  boundary,  and  Impingement 
shock.  In  actuality,  this  region  may  not  be  Isentroplc  as  assumed. 

The  expansion  waves  generated  by  the  expansion  of  the  gas  from  the 
nozzle  exit  pressure  to  the  base  pressure  reflect  from  the  Jet  boundary 
as  compression  waves.  These  waves  will  coalesce  to  form  a boundary 
shock  at  some  finite  distance  downstream  of  the  nozzle.  However,  if 
the  Impingement  shock  Intervenes  before  a strong  boundary  shock  can 
develop,  the  Isentroplc  assumption  should  be  valid.  This  was  found  to 
be  valid  for  the  cases  considered  in  the  theoretical  investigation. 

Once  the  Jet  boundary  Intersects  the  tube  wall,  an  impinge- 
ment shock  develops  as  a result  of  the  turning  of  the  gases  along  the 
wall.  The  flow  aft  of  the  Impingosent  shock  (the  axially  symmettic, 
rotational  flow  region)  is  highly  rotational  and  must  be  treated  as 
nonlsentroplc.  From  the  launcher  design  viewpoint,  this  is  the  most 
Important  region  because  maximum  wall  pressure  exists  behind  the 
Impingement  shock. 
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Governing  Equations 


The  flow  model  is  amenable  to  existing  f.nalytical  techniques 
which  must  be  solved  by  numerical  methods.  Three  separate  analyses, 
one  for  each  region,  are  used.  The  conical  flow  region,  the  simplest 
to  analyze,  is  governed  by  the  Isentropic  relations  and  geometry  with 
regard  to  the  stream  tube  relation: 
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The  other  two  regions  are  treated  by  the  method  of  characteristics  for 
the  descr  '.bed  flows.  The  compatibility  relations  for  axisymmetric 
flows  involves  the  Independent  position  variable,  r.  They  cannot  be 
solved  in  closed  form  as  is  the  case  for  two-dimensional  plane  flows, 
but  can  be  solved  by  a numerical  technique  where  the  terms  which 
involve  the  variable,  r,  are  considered  as  correction  terms  which 
modify  the  two-dimensional  solution  and  adjust  it  to  the  axially 
symmetric  case. 

The  impingement  shock  emanating  from  the  wall  separates  the 
two  axially  sjr  -letric  regions.  The  well  known  oblique  shock  wave 
equations  are  used  in  conjunction  with  the  characteristics  equations 
to  define  the  flows  across  the  shock  and  to  define  the  shape  and 
position  of  the  shock.  An  added  complexity  exists  in  the  solution  of 
the  compatibility  equations  for  the  rotational  flow  region  because  the 
shock  is  curved  and  the  entropy  downstream  of  the  shock  varies  from  one 
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streamline  to  another.  The  entropy  variation  oust  be  taken  ln*-o  account 
In  the  numerical  technique  used  to  generate  the  theoretical  flow  field. 

The  finite  difference  equations  for  steady,  axlsymmetrlc, 
potential  flows  and  for  steady,  axlsymmetrlc,  rotational  flows  are 
presented  In  Appendix  A.  These  equations  are  based  upon  the  method 
characteristics  for  these  flows  and  are  readily  solved  on  computers. 

Computation  Scheme 

A computer  program,  based  upon  this  computation  scheme  to 
solve  the  governing  equations  for  the  flow  model  of  Figure  5,  has  been 
developed.  As  Is  the  nature  of  the  method  of  characteristics,  the  com- 
puter program  provides  a step-by-step  output  which  completely  describes 
the  flow  from  the  c'lnlcal  nozzle  exit  to  the  tube  wall.  Because  of  an 
approximation  which  was  incorporated  Into  the  axially  symuetrlc,  rota- 
tional flow  region,  all  the  solutions  In  this  region  cannot  be  consid- 
ered exact.  The  numerical  results  of  the  analytical  Investigation  were 
obtained  using  the  CDC  6600  computer  system  at  Redstone  Arsenal, 

Alabama. 

The  numerical  routine  Input  Includes  the  ratio  of  specific 
heats  for  the  gas,  V.  and  the  nozzle  geometry,  r*,  r^,  and  6^,  and  cal- 
culates the  flow  conditions  at  the  lip  of  the  nozzle  exit.  The  pressure 
In  the  dead  air  region  or  base  pressure  Is  an  input.  This  pressure  must 
be  obtained  from  actual  tests  or  from  an  ancillary  computation  method 
and  must  be  less  than  the  exit  pressure.  The  Mach  number  la  determined 
using  the  base  pressure/total  pressure  ratio  In  the  isentroplc  pressure 
relation  and  the  turning  angle  of  the  flow  Is  determined  using  the 


19 


Prandtl-Meyer  function.  The  flow  turning  angle  at  the  corner  is 
divided  into  equal  increments,  and  the  values  (v  -f  6)  and  directions 
(6  - ;i)  of  the  right  running  characteristics  which  define  the  expansion 
fan  are  calculated  for  each  increment.  The  flow  angle  and  Mach  number 
of  the  Invlscld  flow  boundary  and  che  right  running  characteristics 
emanating  from  the  nozzle  lip  are  thus  determined. 

Now  a new  position  in  the  conical  flow  region  side  of  the 

expansion  fan  i.s  determined  by  projecting  a radius,  R^,  of  angle, 

0 - A0,  (where  d0  1s  small)  unt..!  it  intersects  the  first  right  run- 

n 

nlng  characteristic.  This  is  a simple  trigometrlc  calculation.  By 
knowing  R^,  the  Mach  number  can  be  calculated  by  the  stream  tube  rela- 
tion and  the  values  (v  - 0)  and  direction  (0  + p)  of  the  left  running 
characteristic  can  be  determined.  This  method  is  also  used  to  deter- 
mine the  characteristics  in  front  of  the  shock  wave,  if  the  shock  wave 
projects  into  the  conical  flow  region  instead  of  the  expansion  fan.  A 
new  right  running  characteristic  is  also  defined  for  this  point.  This 
procedure  is  relatively  simple  and  provides  the  same  solution  as  con- 
structing an  axisymmetrlc  characteristics  net  from  the  nozzle  exit  sur- 
face to  the  expansion  fan  or  shock,  whichever  is  adjacent  to  the 
conical  flow  region. 

The  left  running  characteristic  is  projected  through  the 
expansion  fan  to  the  jet  boundary  by  the  method  of  characteristics  for 
axially  symmetric,  potential  flow.  At  e.''ch  point  in  the  flow  and  on 
the  jet  boundary,  new  right  running  characteristics  are  determined. 

The  computation  method  used  is  taken  from  Owczarek  [24].  The  left 
running  characteristics  are  projected  from  the  conical  flow  region  side 
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of  the  expansion  fan  until  the  Jet  boundary  Intersects  the  tube  vail. 
The  Jet  boundary  flow  characteristics  (l.e.,  primarily,  the  flow  angl» 
because  the  Mach  number  Is  determined  from  the  base  pressure)  at  the 
wall  Intersection  are  determined  by  Interpolation.  As  the  Jet  boundary 
reaches  the  vail  It  Is  assumed  to  turn  parallel  to  the  wall,  thus 
creating  an  oblique  shock.  The  oblique  shock  equations  are  solved  to 
determine  the  Mach  number  behind  the  shock  at  the  impingement  point 
and  the  total  pressure  drop  across  the  shock  wave  for  the  wall  stream- 
line (dS  * 0 for  a streamline  on  either  side  of  the  shock  wave).  The 
shock  wave  angle  with  respect  to  the  wall  Is  also  determined. 

A new  left  running  characteristic  is  Initiated  at  the  expan- 
sion wave  and  Is  projected  Into  the  potential  flow  region  until  the 
Impingement  shock  Is  Intercepted.  The  conditions  in  front  of  the  shock 
are  determined  by  Interpolation.  Figure  6 shows  a construction  of  the 
characteristics  net  In  front  of  the  Impingement  shock. 

A new  shock  angle  and  the  conditions  behind  the  shock  and 
along  the  wall  in  the  rotational  flow  region  are  determined  by  the 
method  detailed  In  Appendix  B.  The  flow  angle  in  the  conical  flow 
region  Is  decreased  by  i^9,  left  running  characteristics  are  projected 
to  the  shock,  and  the  conditions  In  the  rotational  flow  region  are  cal- 
culated until  6^  - mA9  ■■  0 (l.e.,  the  nozzle  centerline  is  reached).  If 
the  Impingement  shock  projects  into  the  conical  flow  region,  the  condi- 
tions in  front  of  the  shock  are  determined  from  the  conical  flow  equa- 
tions because  the  streamline  does  not  pass  through  the  potential  flow  ' 
region.  A characteristics  net  In  the  rotational  flow  region  is  shown 
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in  Figure  7.  In  reality > the  characteristics  are  curved,  but  because 
of  the  method  used,  only  the  average  angle  is  shown. 

In  the  rotational  region,  a comp..e"e  characteristics  net  is 
not  used.  Because  of  the  entropy  gradient  normal  to  the  streamlines, 
the  streamlines  must  be  followed  to  fully  evaluate  the  flow  and  to  con- 
struct a complete  characteristics  net.  To  simplify  the  computation  In 
this  region,  the  characteristics  net  from  the  wall  streamline  [dS  * 0, 

6 ■ 0,  V ■ f(z)]  to  the  streamline  at  the  shock  is  constructed.  This 
Involves  a calculation  technique  that  is  simpler  than  following  the 
streamline  and  constructing  a complete  net  but  provides  results  which 
are  not  as  exact.  The  accuracy  of  the  output  degrades  proportionally 
with  distance  along  the  wall  downstream  from  the  shock  or  along  the 
shock  as  the  tube  axis  is  approached.  The  accuracy  is  discussed  in 
detail  in  Appendix  B.  Because  of  this  simplification,  the  flow  condi- 
tions in  the  region  behind  the  shock  cannot  be  determined  to  the  tube 
axis,  and  the  reflected  shock  or  conditions  behind  it  cannot  be 
ascertained.  The  pressure  immediately  oehind  the  impingement  shock  is 
"exact"  according  to  the  numerical  technique  and  assumptions.  Beca-ise 
the  maximum  wall  pressure  is  assumed  to  be  that  Immediately  behind  the 
impingement  shock,  an  accurate  calculation  of  this  pressure  is  necessary 
for  the  objectives  of  this  study. 

It  was  hoped  that  some  Insight  into  the  flow  aft  of  the  shuck 
could  also  be  obtained  through  this  analytical  technique. 
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Rotational  Flow  Region 


CHAPTER  IV 


EXPERIMENTAL  INVESTIGATION 
General 

A total  of  43  cold  flow  and  8 hot  flow  (actual  static  rocket 
firings)  tests  were  conducted  to  study  experimentally  the  supersonic 
flow  field  In  the  cylindrical  tubes.  These  tests  consisted  of  exhaust- 
ing an  onderexpanded , supersonic  Jet  (l.e.,  P > P ) from  a ^tationary 

e a 

conical  nozzle  Into  a 2.55S-lnch  internal  diameter  tube.  A tually,  a 
total  of  seven  tubes  were  used  to  simulate  the  geometry  wlt.i  the 
rocket  at  various  stations  in  the  launcher.  The  tubs  leng  -hs  ranged 
from  0.7  to  6.2  tube  Internal  diameters.  Since  the  assump  :ions  of  the 
theoretical  analysis  more  closely  represent  the  cold  gas  experiments, 
the  majority  of  the  tests  conducted  were  cold  gas  tests. 

The  nominal  diameter  of  the  nozzle  exit  for  each  type  test 
was  2.1  Inches  which  provides  a relatively  small  ratio  of  tube  radius 
to  nozzle  exit  radius,  x^/r^  1.2.  Test  times  were  relatively  short, 

l.e.,  approximately  10  mlllijeconds  for  the  rocket  firings  and  300 
milliseconds  for  the  cold  gas  tests.  High  frequency  response  pressure 
transducers  were  used  to  determine  the  flow  characteristics  within  the 
tube.  Of  particular  Interest  from  the  launcher  design  viewpoint  were 
the  maximum  wall  pressures.  Static  wall  pressures  were  obtained  along 
the  tube  length.  Special  effort  was  made  to  locate  instrumentation  at 
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Che  intersection  of  Che  shocks  with  the  well  and,  thus,  obtain  maximum 
wall  pressures  and  an  insight  Co  Che  shock  structure. 
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Pitot  pressures  were  recorded  during  the  cold  flow  tests  by 
insertion  of  a pitot  rake  into  the  aft  portion  of  Che  tube.  By  varying 
Che  length  of  Che  Cubes,  the  flow  may  be  probed  in  the  longitudinal 
flow  direction.  The  flow  was  also  probed  in  the  transverse  direction. 
The  pitot  rake  was  used  primarily  to  experimentally  characterize  the 
nozzle  flow  in  Che  conical  flow  region.  High  speed  pictures  were 
obtained  of  the  flow  exiting  from  the  Cube  Co  aid  in  the  analysis  of 
Che  flow  within  Che  Cube. 

Test  Apparatus  and  Instrumentation 

A photograph  of  a typical  cold  gas  test  setup  is  shown  in 
Figure  8.  Note  Chat  the  pitot  rake  is  positioned  outside  Che  tube. 

Pitot  pressures  were  not  obtained  during  the  hot  flow  tests.  The  black 
and  white  striped  board  in  Figure  8 is  a background  for  Che  cameras. 

The  camera  coverage  for  both  type  tests  (cold  and  hot  flow)  consisted 
of  two  cameras,  a high  speed  16-mm  Fastax  camera  running  at  3000  Co 
4000  frames/second  and  a 70-mm  Hulcher  camera  sequencing  at  20  frames/ 
second.  Color  film  was  used  in  both  cameras. 

The  cold  gas  facility  with  nozzle  and  tube  in  place  is  shown 
in  Figure  9.  The  facility  consists  essentially  of  the  motor,  double 
burst  diaphragms,  reservoir,  and  chamber.  Before  operation,  the  reser- 
voir and  chamber  volumes  (0.5  cubic  foot)  are  charged  by  remote  control 
with  dry  nitrogen.  Instrumentation  for  the  operational  control  of  the 
facility  consists  of  four  strain-gage  pressure  transducers  (Baldwin 

26 


Figure  9.  Cold  Gas  Facility  with  Nozzle  and  Tube  in  Place 


gages  having  a 5000-pslg  maximum  pressure)  which  sense  pressures  In  the 
reservoir.  In  the  chamber.  In  the  space  between  the  burst  diaphragms, 
and  In  the  motor  section.  Tests  are  Initiated  by  a double-diaphragm 
technique.  The  gas  In  the  reservoir  and  chamber  then  flows  through 
the  nozzle  like  a blow-down  wind  tunnel. 

Initiation  by  the  double-diaphragm  method  was  accomplished  as 

follows: 


The  reservoir  and  chamber  were  pressurized  to  a 
4020  *20  psla.  At  the  same  time,  the  space  between  the  burst  dia- 
phragms Is  pressurized  to  2020  *20  psla  so  that  at  no  time  will  the 
pressure  differential  across  either  diaphragm  be  greater  than  the  dia- 
phragm burst  pressure.  The  gas  between  the  diaphragms  Is  then  evacu- 
ated Into  the  motor  until  the  pressure  differential  across  the  first 
diaphragm  exccedn  the  rupture  pressure.  The  second  diaphragm  then 
ruptures  allowing  blow  down  through  the  nozzle. 

The  scored  stainless  steel  burst  diaphragms  were 
designed  and  obtained  commercially.  The  diaphragm  rupture  pressure 
was  between  2375  and  2625  psl. 

The  conical  convergent-divergent  nozzle  Is  attached  to  the 
short  motor  section  which  Is  located  aft  of  the  burst  diaphragms 
(Figure  9) • Notice  the  "0"  ring  which  Is  used  to  prevent  leakage 
between  the  nozzle  and  the  forward  section  of  the  tube.  In  the  case  of 
the  hot  gas  tests,  cwo  ”0"  rings  were  used.  It  was  noted  by  Goethert 
[8],  that  leakage  around  the  nozzle  affects  the  base  pressure.  The 
rockets  used  In  the  hot  gas  tests  had  a conical  convergent -divergent 
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nozzle  similar  to  the  nozzle  used  in  the  cold  gas  tests.  The  nozzle 
geometry  and  gas  characteristics  for  both  nozzles  are  presented  in 
Table  I. 


Table  I.  Nozzle  Geometry  and 
Gas  Characteristics 


Parameter 

Cold  Gas 

Rocket 

Nozzle 

Nozzle 

7 

1.4 

1.23 

Throat  Diameter 
(inches) 

1.126 

0.865 

Exit  Diameter 
(inches) 

2.084 

2.065 

Divergence  Cone 

half-angle 

(degrees) 

20 

20 

Mach  Number 
Exit 

2.777 

2.941 

0.03815 

0.02493 

A typical  hot  gas  setup  is  depicted  in  Figure  10.  Double- 
base solid  propellant  rocket  motors  were  used  in  the  hot  gas  tests. 
These  particular  motors  burn  efficiently  at  a temperature  of  approxi- 
mately 5000°F.  These  qualities  (efficient  burning  at  a relatively  low 
temperature),  in  conjunction  with  the  short  burning  time,  minimize  the 
heat  transfer  and  real  gas  effects  (which  include  afterburning). 

The  motor  head  closure  was  replaced  for  the  static  tests  by  a 
pressure/thrust  adapter  shown  in  Figure  10.  The  adapter  facilitated 
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static  firing  and  the  measurement  of  the  motor  chamber  pressure  which 
was  sensed  by  a 5000-pslg  maximum  pressure  Baldwin  gage.  The  existing 
nylon  nozzle-closure  was  replaced  with  an  friable  closure  constructed 
of  styrofoam  to  eliminate  flow  interference  caused  by  the  closure  pass- 
ing through  the  tube. 

Seven  tubes  were  manufactured  according  to  the  drawing  In 
Figure  11.  The  station  dimensions  in  Figure  11  were  non-dlmenslonallzed 
In  terms  of  the  tube  Inside  diameter  and  were  referenced  to  the  nozzle 
exit  plane.  The  ports  along  the  top  of  the  tubes  were  provided  for  the 
wall  static  pressure  transducers.  Section  A Indicates  the  stations 
where  wall  static  pressures  are  taken  around  the  circumference  of  the 
tubes.  A typical  static  pressure  port  Is  also  shown.  A 1/8-lncl 
diameter  by  0. ISO-Inch  long  port  separates  the  transducer  pressure 
sensing  surface  from  the  tube  wall.  The  edges  of  the  orifices  were 
debarred  to  minimize  flow  Interference.  The  relatively  large  port 
diameter  was  needed  to  maintain  the  required  response  time.  Two  types 
of  strain-gage  pressure  transducers,  l.e.,  2000-pslg  maximum  pressure 
Photo-Cons  and  500-pslg  maximum  pressure  Baldwins,  were  used  to  obtain 
tube  wall  static  pressures.  These  gages  were  screwed  Into  the  threaded 
portion  of  the  port.  Sealing  was  accomplished  by  torqulng  the  gage 
bodies  against  an  aluminum  ring.  The  enlarged  v Id  Immediately  in 
front  of  the  pressure  sensing  surface  was  less  than  0.005-lnch  long. 

During  the  course  of  testing.  It  was  found  that  static 
pressure  data  were  required  at  other  positions  along  the  tube  wall. 
Subsequently,  the  tubes  were  modified  to  Incorporate  additional  pres- 
sure transducers.  Table  II  lists  these  additional  pressure  port 
i 
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positions.  CcnsequenCly,  two  series  of  tubes  are  designated.  Series  A 
consists  of  the  tubes  with  port  positions  as  depicted  in  Figure  11,  and 
series  B consists  of  the  same  tubes  with  the  added  ports  according  to 
Table  II. 


Table  II.  Additional  Tube  Wall  Static 
Pressure  Positions 
(Series  B) 


Tube 

z Position 
(ID) 

0 

(deg) 

I 

0.1 

60 

0.3 

120 

3 

1.5 

90 

4 

2.25 

135 

2.50 

45 

2.75 

180 

5 

3.25 

135 

3.50 

45 

3.75 

180 

6 

4.25 

135 

4.50 

45 

4.75 

180 

7 

5.25 

135 

5.50 

45 

5.75 

180 

Two  views  of  the  pitot  rake  in  situ  are  shown  In  Figure  12. 

The  adjustable  mount  allows  the  rake  poaltlon  to  be  varied  relative  to 

I 

the  tube  In  the  axial  and  radial  directions.  Mote  that  in  this  drawing 
the  probes  are  positioned  within  the  tube.  The  rake  consists  of  three 
square  faced  1/8-lnch  outside  diameter  cylindrical  Cubes  mounted  along 
the  edge  of  a 17-1/2-degree  half  angle  wedge.  The  diameter  of  the  hole 
In  the  face  of  each  probe  was  0.040  Inch  and  provides  a ratio  of  hole 
diameter  to  cylinder  outside  diameter  of  0.32.  The  Cubes  are  posi- 
tioned horizontally  as  shown  and  project  0.61-lnch  forward  from  the 
wedge  edge.  A shroud  Is  attacked  to  protect  the  pressure  transducers 
which  screw  Into  Che  aft  end  of  Che  wedge.  The  pitot  pressures  were 
sensed  by  2000-palg  maximum  pressure  Photo-Con  strain-gage  type  pres- 
sure transducers.  The  Cube  distance  Co  Che  pressure  sensing  surface  of 
Che  gage  was  1.61  Inches  with  Che  enlarged  void  Immediately  in  front  of 
the  pressure  sensing  surface  being  less  than  0.005-lnch  long.  These 

distances  were  kept  at  a minimum  to  Insure  adequate  time  response.  . 

Sealing  was  accomplished  by  torqulng  the  gage  bodies  against  an  alumi- 
num ring. 

I 

Test  Description 


Pertinent  Information  regarding  each  test  Is  presented  in 
Table  III.  The  Cube  wall  static  pressure  locations  are  not  given. 
However,  In  every  case,  wall  static  pressures  were  recorded  at  each 
port  location.  These  locations  for  each  Cube  are  defined  in  Figure  11 
and  Table  II.  By  slightly  moving  Che  Cube  in  the  axial  direction,  all 
Che  port  locations  relative  to  Che  nozzle  exit  are  changed  such  as  for 
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Table  III.  Experimental  Program  Summary 


Test 

Cold 

Hot 

Tube 

Tube 

Position 

Pitot  Position 

Preliminary  Tests 

1 

X 

No 

NA 

PNU 

2 

X 

No 

NA 

PNU 

3 

1 

X 

No 

NA 

Nozzle  exit 

4 

X 

No 

NA 

Nozzle  exit 

5 

X 

la* 

Normal 

0.5-ln.  In  tube 

6 

X 

la* 

Normal 

Tube  exit 

7 

X 

la* 

Normal 

PNU 

8 

X 

3a* 

Normal 

Tube  exit  plane 

9 

X 

3a* 
Tes  tt 

Normal 

3 

Tube  exit  plane 

10 

X 

No 

NA 

Nozzle  exit 

11 

X 



No 

NA 

0.22-ln.  aft 
nozzle  exit 

* Not  sealed  between  nozzle  PNU  Pl^ot  not  us>  1 

and  tube  NA  fiot  applicable 
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Table  III.  Continued 


- 

Test 

Cold 

Hot 

Tube 

Tube 

Position 

Pitot  Position 

12 

X 

No 

NA 

1.281-in.  aft 
nozzle  exit 

13 

X 

No 

NA 

2.558-in.  aft 
nozzle  exit 

14 

X 

No 

NA 

5.116-in.  aft 
nozzle  exit 

15 

X 

i 

No 

1 

NA  1 

7. 674-in.  aft 
nozzle  exit 

16 

X 

No 

NA 

:-0.232-ln.  aft 
nozzle  exit 

17 

X 

No 

NA 

12.790-ln.  aft 
nozzle  exit 

18 

X 

No 

NA 

15.348-in.  aft 
nozzle  exit 

19 

X 

1 

1 . 

: 1 
1 ' 

la 

0.25-in. 

aft 

0.5  in.  In  tube 

20 

X 

1 

\ 

la 

Normal 

0.5  in.  In  tube 

21 

X 

1 

1 

la 

i 

Normal 

0.5  in.  in  tube 

22 

X 

la 

Normal 

0.5  in.  in  tube 

f ■ i 

NA  -•  Not  applicable  • 
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Table  III.  Continued 


Test  Cold  Hot  Tube 


1 

Tube 

Tube 

Position 

Pitot  Position 

! 

1 

Normal 

Outside  tube 

7. 

! 

Normal 

0.5  In.  In  tube 

7a 

Normal 

Tube  exit 

' 

Normal 

0.5  In.  in  tube 

2a 

j 

Normal 

0.5  In.  In  tube 

i 

2a 

Normal 

0.5  In.  In  tube 

1 

2a 

Normal 

Outside  tube 

3a 

Normal 

0.5  In.  in  tube 

3a 

Normal 

Outside  tube 

4a 

Normal 

0.5  In.  In  tube 

4a 

Normal 

Outside  tube 

5a 

Normal 

0.5  In.  In  tube 

5a 

I 

1 Normal 

1 

Outside  tube 

6a 

! Normal 

0.5  In.  In  tube 

6a 

Normal 

Outside  *’ube 
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Table  III.  Continued 


Test 

Cold 

Hot 

Tube 

Tests 

Tube 

Position 

Pitot  Position 

38 

X 

3b 

Normal 

PNU 

39 

X 

-4b 

Normal 

PNU 

AO 

X 

4b 

0.263 
In.  aft 

PNU 

A1 

X 

5b 

Normal 

PNU 

42 

X 

6b 

Normal 

PNU 

A3 

X 

6b 

0.263 
in.  aft 

PNU 

AA 

X 

7b 

Normal 

PNU 

45 

X 

lb 

Normal 

PNU 

A6 

X 

lb 

0.406 
In.  aft 

Tube  exit 

47 

X 

No 

NA 

Nozzle  exit 

48 

X 

No 

NA 

1.0-ln.  aft 
nozzle  exit 

49 

X 

lb 

Normal 

PNU 

PNU  — Pitot  not  used 
NA  — Not  applicable 
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Tabic  III.  Concluded 


Test 

Cold 

Hot 

Tube 

Tube 

Position 

Pitot  Position 

50 

X 

lb 

0.25-ln. 

aft 

PNU 

51 

X 

4b 

Normal 

PNU 

52 

X 

7b 

Normal 

PNU 

53** 

X 

lb 

Normal 

PNU 

5/,** 

X 

2b 

Normal 

PNU 

55 

X 

lb 

Normal 

PNU 

56 

M 

lb 

0.125  In. 
aft 

PNU 

**  Only  camera  daCa  obtained 


PNU  Pitot  not  used 
Mot  applicable 


/ 
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test  19.  This  was  done  to  provide  additional  pressure  data  and  to 
define  the  pressure  distribution  without  the  addition  of  more  orifices 
and  Instrumentation.  The  pitot  probe  positions  In  the  flow  direction 
are  given.  However,  the  transverse  positions  of  the  pitot  probes  are 
not  given  but  will  be  Indicated  when  the  data  are  presented. 

Preliminary  tests  were  conducted  to  check  out  the  operation 
of  the  cold  gas  facility  and  Instrumentation.  Note  that  for  these 
tests,  the  "0"  ring  was  not  positioned  between  the  nozzle  and  tube  wall 
to  seal  the  dead  air  region  from  the  atmosphere.  These  tests  also  pro- 
vided Information  regarding  the  magnitude  and  variation  of  static  pres- 
sures along  the  tube  wall.  This  Information  was  used  to  calibrate 
the  pressure  gages  over  the  ranges  Indicated  so  that  the  best  accuracy 
could  be  obtained. 

Two  rockets  were  static  fired  to  determine  If  reliable  cham- 
ber pressure  could  be  obtained.  To  obtain  pressure  at  the  head  end  of 
the  rocket  motor,  the  propellant  pin  plate,  which  Is  between  the  pro- 
pellant and  the  pressure  gage,  had  to  be  filled  with  holes.  Of  concern 
was  the  effect  of  this  restriction  on  obtaining  adequate  pressure  time 
response.  Upon  comparison  of  thrust  data  with  the  pressure  data.  It 
was  found  that  this  modification  was  adequate.  The  friable  nozzle 
closure  was  also  checked  at  this  time. 

A total  of  13  cold  gas  tests  were  conducted  with  the  nozzle 
exhausting  Into  the  atmosphere,  l.e.,  without  tubes.  These  exhaust 
flows  were  probed  with  the  pitot  rake  at  various  positions  In  the  flow. 
This  Information  was  used  to  ascertain  the  type  of  flow  provided  by  the 
conical  nozzle. 
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The  remainder  of  the  testa  were  conducted  with  tubes.  Test 
density  per  type  test  per  tube  Is  presented  in  Table  IV.  Pressure  data 

Table  IV.  Test  Density  Per  Type 
Test  Per  Tube 


1 

1 

Tube 

Hot  Gas 

Cold  Gas 

No.  of 
Teats 

No.  of  Tests 

No 

Pitot 

. ' 1 

No 

Pitot 

Pitot 
in  Tube 

Outside 

Tube 

1 

B 

A 

2 

B 

A 

3 

B 

A 

B 

A 

5 

B 

A 

6 

B 

TOTALS 


NOTE:  Atmospheric  tests;  13  cold  gas,  2 hot  gas 


for  tests  53  and  54  were  not  obtained  as  these  tests  were  conducted  to 
obtain  better  pictorial  coverage  of  the  impingement  shock.  Most  of  the 
tests  were  conducted  with  tube  1.  This  tube,  which  was  the  shortest 
and  for  which  the  pressure  ports  were  concentrated  near  the  jet 
boundary-wall  Intersection,  allows  the  flow  near  the  Impingement  shock 
to  be  investigated  with  the  pitot  rake.  Two  other  pressure  peaks  A'ore 
noted  in  the  cold  gas  flows.  Hence,  numerous  tests  with  tubes  4 and  6 
were  conducted  to  define  the  flow  conditions  or  causes  of  the  peaks. 

To  give  a better  wall  pressure  representation,  the  A-series 
tubes  were  altered  to  provide  more  data  points.  Tests  38  through  56 
were  conducted  with  the  altered  tubes,  designated  B series , The  Hulcher 
sequencing  camera  was  incorporated  into  the  instrumentation  for  these 
tests.  Because  of  the  limitation  of  available  rockets,  only  eight  hot 
gas  tests  were  conducted  including  the  two  preliminary  tests.  Conse- 
quently, most  of  the  effort  was  directed  toward  obtaining  wall  pres- 
sures in  tube  1.  Although  tests  were  run  with  tubes  4 and  7 to  obtain 
an  idea  of  the  pressure  conditions  along  the  wall,  the  test  times  for 
the  hot  gas  tests  were  too  short  to  obtain  70-mm  camera  coverage. 

Data  Reduction  and  Discussion  of  Test  Data 

For  the  cold  flow  tests,  the  pressures  sensed  by  the  trans- 
ducers were  recorded  on  tape  and  on  Consolidated  Electronic  Corporation 
oscillograph  recorders.  The  recorders  were  run  at  160  inches/second. 

The  gages  were  calibrated  in  the  required  pressure  range  before  a 
series  of  tests.  The  maximum  pressure  to  be  sensed  by  each  gage  was 
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estlnated  and  this  presaure  was  establlahed  as  the  required  range. 

This  range  was  projected  as  full  scale  on  the  recorder  to  obtain  the 
maximum  resolution  possible. 

The  data  recorded  on  the  tape  was  transcribed  from  analog  to 
digital  form  so  that  It  could  be  reduced  to  engineering  units  by  the 
CDC  6600  computer.  The  computer  program,  which  was  written  to  reduce 
the  data,  al$o  performed  necessary  calculations  (viz,  averaging,  non-* 
dlmenslonallzlng,  and  Mach  number  calculations) . To  obtain  the  best 
accuracy  possible,  the  zero  shifts  which  were  noted  on  the  oscillograph 
records  were  Inputed  to  the  computer  program  for  correction.  Because 
the  velocities  in  the  cold  gas  facility  reservoir  and  chamber  were 
small,  Che  average  of  these  two  pressures  were  taken  as  the  total  pres- 
sure. All  pressures  were  nondlmenslonalized  by  the  total  pressure. 

The  double  burst  diaphragm  technique,  which  Is  also  used  In 
shock  Cubes,  Initiates  a series  of  shocks  and  expansion  waves  which 
propagate  Into  the  tube.  The  time  lag  of  the  Instruments  when  checked 
against  these  pressure  pulses  was  nil.  It  was  decided  Co  reduce  Che 
data  over  that  portion  of  the  test  which  represented  a near  steady 
state  and  which  was  void  of  pressure  discontinuities.  Consequently, 

Che  data  were  analyzed  over  Che  time  span  of  0.023  to  0.073  second 
after  flow  Initiation.  The  total  pressure  and  typical  wall  pressure 
for  test  19,  showing  time  Interval  between  reduced  points  and  time 
span  for  which  data  were  reduced,  are  presented  In  Figure  13. a.  The 
computer  program  provided  data  over  this  time* span  in  Increments  of 
2.5  milliseconds.  Although  Che  total  test  run  time  was  approximately 


45 


0.3  second,  the  data  were  not  reduced  after  0.073  second  to  eliminate 
the  reduction  of  data  when  condensation  was  present  In  the  flow. 

Based  upon  the  literature  survey.  It  is  reasonable  to  assume 
that  the  conditions  In  the  tube  are  Independent  of  the  nondimens lonal~ 
Izing  parameter,  P^.  With  this  In  mind  and  for  the  ease  of  data  analy- 
sis, only  cold  flow  data  corresponding  to  five  total  pressures  are 
considered,  i.e.,  ■>  2300,  2000,  1700,  1500,  and  1300  psla.  Because 

the  reduced  data  may  not  provide  Inforization  at  these  exact  total  pres- 
sures, the  data  presented  represent  an  average  of  three  total  pressures 
one  on  each  side  and  the  closest  total  pressure  (Figure  13. b).  Note 
the  data  to  be  averaged  are  already  nondlmenSlonalized.  The  time  span 
over  which  these  three  pressures  are  recorded  Is  5 milliseconds  which 
provides  an  adequate  representation  of  the  flow  field  for  this  time. 

The  use  of  the  average  of  the  three  values  inherently  tends  to  smooth 
the  data.  The  Mach  numbers  were  averaged  as  were  the  static  wall  pres- 
sures when  a large  variation  In  Mach  number  was  noted.  However,  in 
most  cases,  the  Mach  number  varied  only  a few  hundredths  and.  In  these 
cases,  the  Mach  number  presented  Is  an  average  of  all  the  values. 

The  use  of  the  tape-computer  system  was  not  feasible  to 
reduce  the  hot  flow  data  because  of  the  short  test  times.  The  hot  flow 
data  were  directly  reduced  from  the  oscillograph  records.  Data  repre- 
senting three  times  thst  occurred  during  the  motor  operation  were 
reduced.  A typical  pressure/tlme  representation  is  shown  in  Figure  14. 
The  pressure  data  were  nondlmenslonalized  with  the  motor  chamber  pres- 
sure. For  hot  gas  experiments  55  and  56,  the  motor  chamber  pressures 
were  markedly  different  from  that  shown  In  Figure  14.  The  pressures 
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AND  A (p^l 


Figure  14.  Typical  Pressure/Time  Representation 


For  Hot  Gas  Test 


peaked  at  approximately  6380  psia  and  were  regressive  In  nature,  over  a 
shorter  time.  This  occurred  because  the  temperature  of  the  propellant 
grains  were  higher  than  those  of  the  other  hot  gas  experiments. 

The  nondlmenslonal  data  represents  two  pressure  measurements  - 
one  divided  by  another  - with  Individual  accuracies  (plus  or  minus  a 
certain  percent)  which  may  be  additive  or  nullifying.  The  accumulative 
error,  or  worst  accuracy.  Is  given.  However,  statistically,  and  in 
reality,  the  average  error  will  be  less.  According  to  Gracey  [25],  a 
pitot  probe  of  this  design,  l.e.,  with  a ratio  of  hole  diameter  to 
probe  outside  diameter  of  0.32,  will  be  insensitive  to  an  angle  of 
attack  of  *13-1/2  degrees.  Because  the  static  pressures  In  the  tube 
vary  as  a function  of  radius,  the  Mach  number  of  the  flow  downstream  of 
the  shocks  cannot  be  determined  from  the  pitot  pressure  and  wall  static 
pressures.  Also,  because  of  design  flow  angle  limitation,  the  probes 
are  limited  to  flow  angles  of  less  than  *13-1/2  degrees.  A knowledge 
of  the  flow  angle  is  required  and  was  not  obtained  experimentally. 

Also,  the  ability  of  the  pitot  probe  to  accurately  measure  pitot  pres- 
sure In  a flow  with  a large  radial  flow  gradient  is  in  question.  The 
overall  accuracy  of  the  pitot  probe.  Including  response  time  errors. 
Instrument  and  positioning  error,  and  the  method  of  data  reduction, 
but  excluding  flow  angle  and  radial  gradient  sensitivities,  is  con- 
sidered to  be  within  *6  percent.  This  amounts  to  a Mach  number  error 


of  approximately  *3  percent.  The  error  for  the  wall  static  pressures 
obtained  In  the  cold  flow  tests  Is  believed  to  be  *3-1/2  percent.  The 
accuracy  of  the  data  In  the  hot  gas  Is  degraded,  primarlxy  as  a result 


of  the  severity  of  the  pressure  and  temperature  gradients  with  respect 
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to  time  imposed  on  the  transducers  and  the  resolution  of  data  during 
reduction.  This  error  can  be  as  high  as  *5  percent.  The  static  pres- 
sure accuracies  are  based  upon  the  response  time,  error  of  the  data 
acquisition  system,  instrument  error,  zero  shift,  errors  Imposed  by 
orifice  edge  form  and  port  diameter,  and  errors  due  to  the  method  of 
data  reduction. 
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CHAPTER  V 


RESULTS  AND  DISCUSSION 
Preliminary  Remarks 

Before  any  conclusions  can  be  reached  concerning  the  nature, 
ot  the  fJow  in  the  tube  and  the  parameters  which  govern  the  flow, 
certain  basic  premises  must  be  investigated  and  verified  to  ascertain 
the  validity  of  the  experimental  methods  and  assumptions  used  in  the 
theoretical  analysis.  Because  of  the  method  of  testing,  l.e.,  multiple 
tests  with  limited  instrumentation  per  test,  the  repeatability  of  the 
flow  conditions  within  the  tube  from  experiment  to  experiment  is 
paraoxDunt . 

The  Initiation  conditions  and  the  time  Interval  over  which 
the  cold  flow  data  were  obtained  were  standardized  in  an  effort  to 
reproduce  the  flows  in  the  Cube  from  experiment  to  experiment.  This 
was  also  attempted  for  the  hoc  gas  rocket  firings;  however,  two  tests 
(55  and  56  as  discussed  previously)  were  inadvertently  conducted  at  a 
different  propellant  conditioning  temperature.  The  total  pressures  for 
these  two  tests  were  markedly  different  and  required  that  the  data  be 
obtained  over  a different  time  span. 

The  correlation  of  the  experimental  results  with  those 
obtained  by  theoretical  analysis  depends  upon  how  well  the  theoretical 
assumptions  were  duplicated  during  the  experimental  Investigations.  The 
character  of  the  nozzle  exhaust  (underexpanded,  supersonic,  conical, 
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and  axlsymmetrlc)  and  the  symmetry  of  the  flow  in  the  tube  were  depen- 
dent on  the  constraints  of  the  experiment  and  are,  to  a certain  extent, 
controllable  by  the  experimenter.  The  validity  of  the  fully  developed 
supersonic  flow  assumption  upon  which  the  flow  model  was  developed  and 
the  quasi-steady  flow  assumption  Is  dependent  upon  the  results  of  the 
experiments. 

Character  of  the  Nozzle  Exhaust 

The  Mach  number  In  the  nozzle  exhaust  was  calculated  from  the 
total  pressures  behind  a normal  shock,  as  measured  by  the  pitot  probes, 
and  the  total  pressure  In  front  of  the  normal  shock,  as  measured  In  the 
reservoir.  The  Mach  numbers  in  the  proximity  of  the  nozzle  exit  were 
used  to  characterize  the  nozzle  flow  and  to  ascertain  If  the  nozzle  did. 
In  fact,  provide  conical  flow.  Because  pitot  data  were  not  obtained, 
during  the  hot  gas  tests,  only  the  cold  gas  nozzle  exhaust  character- 
istics can  be  presented. 

Figure  15  presents  the  Mach  number  obtained  from  the  pitot 
probes  plotted  with  respect  to  the  physical  conical  radius  parameter, 
(R-Zg)/r^.  Zero  represents  the  point  on  the  nozzle  centerline  In  the 
exit  plane.  Data  are  presented  for  the  nozzle  exhausting  into  the 
atmosphere  and  into  tubes  1 and  2,  in  the  vicinity  of  the  nozzle  exit. 
Also  Included  is  the  theoretical  curve  for  conical  flow.  Note  that 
the  data  points,  with  the  exception  of  the  numbered  points,  fall  along 
the  theoretical  curve.  The  physical  positions  of  the  probes  are  shovm 
In  Figure  16.  Figure  16. a shows  the  salient  features  obtained  from 
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the  theoretical  analysis  (P]j/Pq  *■  0.0218)  which  define  the  conical 
flow  region  for  the  nozzle  exhausting  In  the  tubes.  The  Inner  edge 
of  the  expansion  fan  which  u<;flnes  the  conical  flow  region  for  the 
nozzle  exhausting  into  the  atcicay*  ere  is  shown  in  Figure  16. b.  The 
numbered  points  are  either  located  on  or  downstream  of  the  Impingement 
shock  outside  of  the  conical  flow  region.  Point  C represents  a pitot 
probe  position  for  test  23.  Figure  17  is  a photograph  taken  from  16-mm 
film  of  this  test  condition.  The  impingement  shock  emanating  from  the 
tube  and  the  normal  shock  Induced  by  the  presence  of  the  pitot  probes 
are  evident  In  this  photograph.  The  Impingement  shock  Is  located 
Inboard  of  the  outboard  pitot  probe,  data  point  C,  and  outside  of  the 
next  probe  as  shown  In  Figure  16. a.  The  Mach  number  data  from  the 
points  outside  of  the  conical  flow  region  are  Invalid  because  of  the 
presence  of  the  shock.  However,  they  serve  to  give  an  indication  of 
the  physical  bounds  to  the  conical  flow  region. 

Thus,  it  is  seen  that  the  nozzle  provides  conical  flow  in  the 
vicinity  of  the  exit.  Because  the  pitot  probes  were  located  at  dif- 
ferent radial  positions  In  Che  nozzle  exhaust,  it  is  also  concluded 
Chat  Che  nozzle  provides  flow  which  Is  a}d.ally  symmetric.  Because  the 
Mach  numbers  as  determined  from  the  probe  data  were  essentially  con- 
stant over  the  entire  data  acquisition  time  span  or  total  pressure 
range,  the  nozzle  exhaust  flow  Is  independent  of  total  pressure  or  time. 

The  constraints  of  the  experimental  control  conditions,  i.e., 

$ 

nozzle  geometry  coupled  with  total  pressures,  dictate  that  the  nozzle 
exhaust  is  underexpanded  with  respect  to  the  tube  exit  pressure,  P^. 

This  does  not  assure  that  the  nozzle  exhaust  into  the  tube  with  the 
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gure  17.  View  Showing  Impingement  Shock  Emanating  from  Tube  and 
Normal  Shock  Induced  by  Pitot  Probes  for  Test  23 


unknown  base  pressure  is  underexpanded.  This  Is  the  case,  however, 
because  the  nozzle  exit  pressure  is  greater  than  the  base  pressure. 
The  base  pressure/nozzle  exit  pressure  relations  obtained  during  the 
cold  and  hot  gas  tests  are  presented  in  Table  V. 


Table  V.  Base  Pressure/Nozzle  Exit  Pressure  Ratios 


Cold  Gas  Tests 

Hot  Gas  Tests 

P 

P. 

P, 

P 

P. 

P 

o 

b 

b 

0 

b 

b 

(Psla) 

' P 

0 

P 

e 

(Psla) 

P 

o 

P 

e 

2300 

.0218 

.571 

5440 

.0066 

.265 

2000 

.0225 

.590 

5350 

.0073 

.293 

1700 

.0231 

.606 

5135 

.0082 

.329 

1500 

.0239 

.626 

1300 

.0244 

.640 

Quasi-Steady  Flow 

According  to  Korst  et  al.  119],  the  nozzle  exhaust  flow  in 
the  viscous  shear  layers  that  bound  the  dead  air  region  along  with  the 
physical  geometric  contralnts  dictate  the  equilibrium  base  pressure  by 
mass  transfer  into  or  out  of  the  dead  air  region.  And  according  to 
Ihrlg  and  Korst  [30],  the  ratio  of  adjustment  time  for  the  transfer  of 
mass  as  compared  to  pressure  disturbance  adjustment  time  in  the  super- 
sonic flow  is  approximately  30  for  Mach  1 flow;  l.e.,  the  pressure 
field  in  the  main  flow  adjusts  almost  instantaneously  to  the  slower 
base  pressure  variation.  Apcordlng  to  Fabrl  and  Slestrunck,  the  ratio 
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Pj^/Pq  for  supersonic  fully  developed  steady  flow  is  a constant.  Korst 
also  predicts  a constant  for  steady  flow,  but  allows  for  a variatlc  in 
his  quasi-steady  analysis. 

It  is  significant  that  the  dimensionless  base  pressure 
Increased  with  time  for  the  cold  and  hot  gas  flows.  However,  there 
were  two  Important  parameters  varlnp  .uce  an  Increase  in  time  was 
accompanied  by  a decrease  in  total  pressure.  Because  the  total  pres- 
sures occurred  at  approximately  the  same  times  during  the  tests  for  the 
cold  flow  data  presented,  the  effects  of  the  two  parameters  could  not 
be  separated.  The  total  pressure  histories  were  also  duplicated  for 
the  majority  of  the  hot  gas  tests.  However,  this  was  not  the  case  for 
two  hot  gas  tests  (tests  55  and  56)  which  were  conducted  at  a higher 
propellant  temperature.  The  burning  rate  of  the  particular  propellant 
varies  with  temperature;  consequently,  the  total  pressures  for  these 
tests  were  higher  and  the  burning  time  was  decreased.  The  base  pres- 
sure is  presented  as  a function  of  total  pressure  for  the  hot  gas  flows 
in  Figure  18.  The  total  pressures  of  the  ordinate  are  reversed  to  give 
an  indication  of  Increase  in  time  from  left  to  right.  This  is  because 
of  the  regressive  burning  characteristics  of  t'  e rocket  motors.  Note 
that  because  of  the  shorter  burning  time  of  tests  55  and  56,  the  time 
derivative  of  the  dimensionless  base  pressure  of  the  two  traces  in 
Figure  18  would  be  approximately  the  same  if  the  dimensionless  back 
pressure  was  considered  with  respect  to  time. 

The  missile  oriented  base  pressure  computer  program,  as 
described  by  Addy  [31],  was  altered  to  calculate  the  base  pressure  for 
the  sudden  expansion  of  an  axisymnetric  jet  into  a shroud  for  the  zero 
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flow  base  bleed  case  by  Walker  [32],  The  calculaclon  mechod  is  based 
upon  Korst's  flow  model.  The  steady-state  base  pressures  for  the  cold 
and  hot  gas  flows  as  obtained  from  [32]  are  presented  in  Table  VI  along 


Table  VI.  Experimental  and  Predicted  Base  Pressure 


Cold  Gas 

Hot  Gas 

P 

P. 

P 

o 

b 

b 

o 

b 

b 

(Psia) 

P 

o 

P 

a 

(Psia)' 

P 

o 

P 

a 

Steady-Stc 

ite  Predict 

ed  Base  Pre 

ssures  [32] 

NA 

1.0244 1 

NA 

NA 

.118 

NA 

Ex' 

>erlmental 

Base  Pressu 

res 

2300 

.0218 

3.41 

5440 

.0066 

2.44 

2000 

.0225 

3.06 

5350 

.0073 

2.66 

1700 

.0231 

2.67 

5135 

.0082 

2.86 

1500 

.0239 

2.44 

1300 

1 l6 

NA  - Not  Applicable 


with  those  obtained  from  experiment.  Also  Included  is  the  base 
pressure/back  pressure  ratio,  which  in  this  case  is  atmospheric 
pressure.  Of  significance  is  the  observation  that  originally  the  base 
pressure  was  equal  to  the  back  pressure  or  atmospheric  pressure  and  is 
increasing  toward  the  equilibrium  value.  Mass  is  being  transfe’rred 
from  the  exhaust  flow  into  the  dead  air  region.  In  the  case  of  the 
cold  flow  tests  where  the  data  are  obtained  over  a much  longer  time 


span  and  where  the  total  preaswre  Is  not  varying  as  markedly  as  In  the 
hot  flow  tests,  the  experimental  base  pressure  does  reach  the  predicted 
equilibrium  value  [32].  It  is  concluded  that  the  base  pressure  depends 
on  the  total  pressure  history  and  that  the  resultant  flow  is  dependent 
upon  the  base  prsssure/total  pressure  ratio.  This  conclusion  is  only 
relevant  for  the  cold  gas  flows  over  the  total  pressure  range  which  are 
void  of  the  unsteady  pressure  waves  and  for  the  hot  gas  flows  over  the 
rocket  motor  operating  chamber  pressures.  The  flows  within  the  tube  can 
be  treated  as  steady  with  the  base  pressure  being  a function  of  time. 

If  the  total  pressure  history  is  known,  the  base  pressure  variation  may 
be  determined  from  a quasi-steady  method  similar  to  reference  [30]. 

Flow  Symmetry,  Flow  Reproducibility,  and  Tube  Length  Effects 

To  obtain  a measure  of  flow  symmetry,  three  static  pressure 
orifices  were  positioned  120-degrees  apart  around  the  circumference  of 
each  tube.  The  symmetry  is  determined  by  an  evaluation  of  the  unifor- 
mity of  the  static  pressures  sensed  around  the  tube.  The  locations  of 
these  orifices  in  each  tube  are  shown  in  Figure  11.  For  tubes  3 through 
7,  the  axial  location  was  the  same.  The  percent  variation  for  these 
three  wall  static  pressures  for  a particular  total  pressure  and  for  a 
particular  experiment  was  determined  by  subtracting  the  smallest  value 
from  the  largest  and  dividing  by  the  average  of  the  three  values.  For 
the  cold  and  hot  gas  tests,  this  was  done  at  each  of  the  five  total 
pressures.  The  maximum  and  minimum  variations  were  noted  and  the  aver- 
age was  calculated  for  each  tube.  These  variations  are  presented  in 
Table  VII. 


61 


Table  VII.  Determination  of  the  Symmetry  of  the  Flow 
Using  Circumferential  Pressure  Distribution 
over  the  Range  of  Stagnation  Pressure 


Location 
from  Nozzle 

A - %\ 

Tube 

Wall  Static  Pressure  Variation 

Maximum 

(percent) 

Average 

(percent) 

Minimum 

(percent) 

1 'e  ) 

Cold  Gas 

0.22 

1 

3.6 

2.9 

1.7 

2.42 

2 

6.4 

4.1 

2.1 

4.82 

3 

6.7 

4.5 

2.1 

4.82 

4 

5.9 

3.2 

0.5 

4.82 

5 

9.7 

5.5 

1.1 

4.82 

6 

6.3 

4.5 

0.1 

Hot  Gas 

1.2* 

1 

3.1 

1.6 

0.4 

1.2** 

1 

7.7 

5.6 

3.1 

4.76 

4 

3.9 

2.4 

0.8 

4.78 

7 

11.1 

8.0 

4.6 

^Nominal  position. 

**These  variations  are  for  high  temperature  rocket  experiments. 


i 

/ 


The  Inherent  error  of  the  static  pressure  measurement  system 
for  the  cold  gas  experiments  Is  *3-1/2  percent.  This  error  represents 
a maximum  variation  of  7 percent.  All  the  variations  for  the  cold  gas 
flows  except  for  tube  5 are  within  this  range.  For  the  hot  gas  experi- 
ments, the  variations  found  for  tube  7 exceed  the  predicted  variation  of 
10  percent  due  to  the  measurement  system.  Under  these  conditions,  It  Is 
difficult  to  distinguish  the  measurement  error  from  an  Indication  of 

unsymmetrlc  flow.  A marked  variation  In  the  static  pressures  measured 
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at  the  three  circumferential  locations  was  noted  at  total  pressures  on 
the  order  of  300  to  400  psia  or  less  in  the  cold  gas  experiments.  The 
indicated  total  pressure  range  provided  nozzle  exit  pressures  which 
were  approximately  equal  to  the  atmospheric  pressure.  It  Is,  therefore 
reasoned  that  these  conditions  are  those  described  by  Fabrl  and 
Selstrunck,  corresponding  to  the  mixed  flow  regime  where  the  Jet 
attaches  to  one  portion  cf  the  tube  wall. 

Of  particular  Interest  In  determining  the  reproducibility  of 
the  flow  In  the  tubes  are  the  variations  of  the  wall  static  pressures 
from  experiment  to  experiment.  These  variations  are  presented  in 
Table  VIII. a along  with  Che  average  static  pressure/coCal  pressure 
ratios  for  the  circumferential  locations.  Because  of  the  small  number 
of  tests  conducted  for  Che  hot  gas  experiments  such  variations  would 
not  have  significance  and  only  values  are  shown.  The  fact  that  the 
variations  are  relatively  small  for  each  tube  Indicates  Chat  the  flow 
was  reproducible  in  the  cold  gas  experiments. 

Also  Included  In  Table  VIII. b are  Che  effects  of  Cube  length 
on  the  static  wall  pressures  at  Che  circumferential  location.  In  Che 
case  of  the  cold  flow  experiments,  the  average  value  and  variations  of 
static  pressures  for  tubes  with  the  same  location  are  provided,  l.e., 
pressures  for  Cubes  3 through  6 with  the  sensor  located  at  the  axial 
position  (z  - Zg)/rg  = 4.82.  The  large  variation  is  a result  of  the 
static  pressures  recorded  for  tube  5 as  evidenced  by  the  variation  with 
this  static  pressure  excluded,  l.e.,  tubes  3,  4,  and  6.  The  reason  for 
the  variation  caused  by  tube  5 is  unexplained.  The  variations  for  the 


Table  Vlll-a.  Variations  of  Static  Pressure  From  Test  to  Teat 
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Total  Pressure  (pal) 

Nominal  values 

Values  are  Cor  'Igh  temperature  rocket  experiments 


hoc  flow  experiments  are  acceptable  with  the  exception  of  the  test^ 
for  which  the  rocket  motor  grains  were  at  the  elevated  temperature 
prior  to  ignition. 

Because  of  the  limited  number  of  pressure  sensors  used  for 
each  Cube,  a clear  representation  of  wall  static  pressure  distribution 
could  not  be  obtained  by  plotting  the  pressure  data  versus  poclclons 
for  each  Individual  tube.  The  data  would  be  too  sparse,  and  conse- 
quently the  slope  of  the  data  curve  at  the  points  would  not  be  deter- 
minable. To  fill  in  the  gaps  and  thus  remedy  this  situation,  the  data 
obtained  for  each  tube  has  Co  be  superimposed  on  Che  same  plot.  This 
can  be  done  only  if  the  flow  is  not  affected  by  the  length  of  t!.e  tube, 
l.e.,  Independent  of  tube  length.  In  the  case  of  f'lly  developed 
supersonic  flow  in  relatively  short  tubes,  the  flow  should  be  indepen- 
dent of  downstream  effects,  Including  tube  length  variations.  The 
only  available  path  for  the  downstream  conditions  to  be  transmitted 
upstream  in  this  type  of  flow  is  through  the  subsonic  portion  of  the 
tube  wall  boundary  layer. 

Because  it  is  not  feasible  to  plot  and  to  identify  every 
pressure  measurement  on  one  ^ aph,  the  nondlmenslonallzed  static 
pressure  data  obtained  at  a particular  location  were  averaged  over  all 
runs  for  the  tube.  In  the  case  of  the  cold  g.  ' experiments,  this  was 
done  for  total  pressures  of  2300,  1700,  and  1300  psla.  These  pressure 
distributions  are  presented  in  Figure  19.  Th^  fact  that  dimensionless 
wall  static  pressures  at  a given  percentage  of  the  cube  length  with 
different  tubes  are  about  the  same  verifies  that  the  flow  is  indepeu- 
dent  of  tube  length.  This  is  also  the  case  for  dimensionless  wall 
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static  pressures  obtained  In  the  hot  gas  experiments.  These  distribu- 
tions for  nominal  total  pressures  of  5440,  5350,  and  5135  psla,  are 
shown  In  Figure  20.  These  nominal  total  pressures  represent  the  range 
of  motor  chamber  operating  pressures  whereas  the  plus  or  minus  values 
shown  In  Figure  20  represent  the  range  for  which  the  static  pressures 
were  determined.  Because  of  the  short  time  span,  the  averaging  method 
which  was  used  for  the  cold  gas  data  was  not  used  for  the  hot  gas  data. 
The  lack  of  data  olnts  In  this  figure  and  the  Incompleteness  of  the 
curve  are  a result  of  the  small  number  of  hot  gas  experiments  which 
were  conducted.  The  observation  that  the  flow  Is  Independent  of  tube 
length  Is  In  consonance  with  that  of  German  et  al.  [15]  which  was  also 
empirically  determined  and  that  of  Fabrl  and  Slestrunck  [3]  which  was 
ascertained  by  theory.  The  reproducibility  of  flow  from  test  to  test 
Is  also  Indicated  by  the  relatively  small  scatter  In  the  data  of 
Figures  19  ard  20. 


Wall  Static  Pressure  Distribution 
The  flow  exiting  the  tubes  was  fully  developed  supersonic 
In  nature  as  evidenced  by  the  shock  wave  and  expansion  waves  which  were 
created  by  the  presence  of  the  pitot  rake  In  the  flow.  This  Information 
was  visible  on  the  16-mm  and  70-mm  films.  Thus  It  Is  concluded  that 
the  flow  was  supersonic  throughout  the  tube. 

The  analytical  flow  model  with  Its  Inherent  assumptions  was 
closely  approximated  In  the  experimental.  Investigation.  To  obtain  a 
correlation  of  theory  with  the  results  of  the  experimental  Investiga- 
tion, the  experimental  base  pressures  together  with  the  experimental 
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geometric  constraints  and  the  ratio  of  specific  heats  were  used  as 
inputs  to  the  theoretical  computer  pregram  for  the  theoretical  analysis. 
This  was  done  for  the  hot  and  the  cold  flows.  The  predicted  wall  static 
pressure  distributions  are  superimposed  on  the  experimental  data  in 
Figures  19  and  20.  As  can  be  seen,  the  theoretical  analysis  adequately 
defines  the  first  pressure  peak.  The  predicted  pressure  data  are  satis- 
factory, indicating  the  approximation  method  used  in  the  theoretical 
analysis  is  satisfactory  with  respect  to  defining  the  pressure 
variations.  The  approximate  theory  developed  herein  cannot  define  the 
second  or  third  peaks.  The  nozzle  exhaust  expands  to  the  base  pressure 
and  the  Jet  boundary-wall  intersection  position  is  a function  of  this 
pressure.  The  pressures  along  the  wall  behind  the  impingement  shock 
are  governed  by  the  base  pressure,  since  it  dictates  the  shock  strength, 
in  essence,  by  dictating  the  jet-boundary  Mach  nuirber  and  flow  angle. 

For  fully  developed  supersonic  flow  in  the  tube,  it  is  reasonable  to 
expect  that  a change  in  base  pressure  will  cause  the  pressures  along 
the  wall  or  the  flow  pattern  in  the  tube  to  change.  However,  note  from 
Figures  19  and  20  that  there  is  a distinct  variation  in  base  pressure 
with  respect  to  total  pressure,  whereas  the  static  pressure  distribu- 
tion along  the  wall  appears  to  change  only  slightly  with  a total 
pressure  change.  As  expected,  the  maximum  wall  static  pressure  occurs 
just  downstream  of  the  impingement  shock-wall  intersection  for  both 
the  cold  and  the  hot  gas  flows  (Figures  19  and  20).  The  maximum  wall 
static  pressures  obtained  from  theory  agree  with  those  found  by 
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experiment.  The  theoretical  maximum  wall  pressure  for  each  base 
pressure  obtained  by  experiment  and  obtained  from  reference  [32]  is 
listed  in  Table  IX. 


Table  IX.  Maximum  Predicted  Wall  Static  Pressure 


Cold  Gas 

Hot  Gas 

**0 

P 

vmax 

PO 

P 

V max 

(Psia) 

^0 

"o 

(Psia) 

^0 

"o 

2300 

0.0218 

0.0870 

5440 

0.0066 

0.0517 

2000 

0.0225 

0.0865 

5350 

0.0073 

0.0526 

1700 

0.0231 

0.0861 

5135 

0.0082 

0.0535 

1500 

0.0239 

0.0854 

1300 

0.0244* 

0.0850 

0.0118  * 

0.0542 

* Predicted  base  pressure  [32] 

1 

Because  the  maximum  wall  pressure  varies  slightly  with  base  pressure, 
an  adequate  estimate  of  this  pressure  can  be  obtained  for  engineering 
purposes  by  using  ths  steady-state  base  pressure  as  predicted  by 
reference  [32] . At  no  time  during  the  experimental  runs  did  the  abso- 
lute value  of  the  maximum  wall  pressure  exceed  that  obtained  at  the 
maximum  total  pressure  during  the  run.  Thus,  it  can  be  concluded  that 
maximum  wall  pressure  occurs  at  the  same  time  as  the  maximum  total 
pressure. 

Supersonic  Flow  Pattern  in  the  Tube 
A further  characterization  of  the  flow  within  the  tube 
hinges  on  the  correlation  between  the  wall  static  pressures  and  the 
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photographic  evidence.  According  to  Fabrl  (Figure  2.c),  the  three 
static  pressure  peaks,  which  are  evident  In  Figures  19  and  20,  are 
caused  by  oblique  shocks  intersecting  the  wall.  The  shock  pattern 
within  the  tube,  as  deduced  from  these  pressure  peaks,  and  the  shock 
locations,  as  manifested  in  the  70-mm  photographs  of  the  flow  exiting 
the  tube,  are  shown  in  Figure  21.  The  impingement  shock  which  causes 
the  first  pressure  peak  is  evident  in  the  70-mm  photographs  of  the 
flow  exiting  from  tubes  1 and  2 (Figures  22  and  23). 

Hie  Impingement  shock  causes  the  flow  in  the  jet  boundary 
to  turn  and  proceed  along  the  tube  wall.  This  shock  emanates  from 
the  jet  boundary-wall  Intersection  and  projects  to  the  tube  axis 
as  an  apparent  paraboloid  surface.  The  time  and,  consequently,  the 
total  pressure  conditions  for  the  flows  shown  in  these  photographs 
are  unknown.  They  are  not  shadowgraphs  or  schlieren,  although  a 
high  intensity  lamp  was  used  for  lighting.  The  process  which  makes 
the  flow  visible  is  the  scattering  of  light  from  the  vapor  condensed 
in  the  flow  itself  as  it  expands  or  on  the  flow  boundary,  in  which 
case,  the  condensed  vapor  was  entrained  from  the  atmosphere.  The 
shock  processes  the  flow  (i.e.,  increases  the  temperature  and 
pressure) , thus  vaporizing  the  condensation  so  that  the  shock  location 
can  be  identified.  This  photographic  technique  also  can  picture  large 
density  gradients  (such  as  a shadowgraph)  as  evidenced  in  some  photos 
which  show  the  shocks  that  develop  because  of  the  pitot  probes.  The 
impingement  shock  Intersects  the  nozzle  axis  just  aft  of  tube  2.  This 
intersection  point  is  indicated  on  the  ordinate  in  Figure  21. 
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Also  shown  In  this  figure  is  the  impingement  shock-axis  location  as 
predicted  by  the  theoretical  analysis.  The  theory  gives  a good 
indication  of  this  point.  However,  it  should  be  noted,  that  in  order 
to  obtain  this  point  by  the  described  theoretical  method,  the  shock 
has  to  be  projected  theoretically  to  the  axis.  The  theoretical 
approach  is  not  expected  '.o  be  accurate  at  the  axis.  Descriptions  of 
the  flow  field  as  obtained  from  the  analytical  investigation  of  the 
cold  and  hot  gas  flows  are  presented  in  Appendix  C.  The  second 
pressure  peak  is  attributed  to  the  impingement  shock  reflecting  from 
the  axis  back  to  the  wall.  This  is  not  apparent  in  the  photograph 
of  the  flow  exiting  tube  3. 

Figure  24  shows  the  shock  projecting  from  the  point  on 
the  wall,  which  corresponds  to  the  second  pressure  peak,  to  the  axis. 
This  photograph  was  taken  at  the  exit  of  tube  4,  The  tube  was 
moved  aft  0.263  inch  in  Figure  24. b.  The  shock-axis  intersection  point 
is  indicated  on  Figure  21.  Both  shock-axis  intersection  points  (im- 
pingement shock  and  second  shock)  occur  when  the  tube  wall  static 
pressures  were  decreasing  to  their  minimums. 

The  exit  planes  of  tubes  3,  5,  and  7 are  near  the  three 
valleys  for  the  wall  static  pressure  and  should  give  an  indication  of 
the  shocks  reflecting  from  the  tube  axis.  No  shocks  were  evident; 
however,  a bulge  or  n expansion  in  the  tube  exhaust  flow  did  occur 
downstream  of  the  tube  exit.  As  it  exits  from  the  tube,  the  flow 
expands  very  little  denoting  that  the  tube  exit  pressure  is  slightly 
larger  than  atmospheric.  At  a finite  distance  downstream,  the  flow 
expands  radically  indicating  a considerable  increase  in  pressure 


which  could  be  caused  by  Che  shock.  The  point  outside  of  Cube  / In 
Figure  21  was  obtained  from  this  bulge.  The  shock  emanating  from  the 
wall  at  the  third  pressure  peak  Is  not  evident  In  r.!  e photograph  of 
the  flow  exiting  Cube  6. 

The  cold  flow  is  supersonic  throughout  the  tubes.  The  flow 
patterns  give  further  evidence  Chat  the  flow  Is  axially  symmetric. 

The  correlation  of  the  pitot  data,  which  were  obtained  at  different 
radial  locations  In  Che  flow,  with  theory  In  the  conical  region  and 
the  analysis  of  Che  circumferential  wall  static  pressures  also  verifies 
that  the  flow  Is  axially  symmetric.  The  anomaly  found  In  Che  evalua- 
tion of  Che  circumferential  static  pressures  for  tube  5 has  not  been 
resolved  and  Is  attributed  to  Instrument  error.  The  16-mm  movie 
coverage  of  Che  hot  gas  experiments  also  revealed  that  the  flows  within 
the  tube  were  supersonic  and  symmetric. 


CHAPTER  V 


CONCLUDING  REMARKS 

The  present  investigation  was  undertaken  to  develop  a design 
model  for  the  flow  field  which  results  when  an  underexpanded  rocket 
exhaust  exits  from  a conical  nozzle  into  a cylindrical  tube.  The 
experimental  phase  of  the  Investigation  employed  short  duration  cold 
gas  flows  and  actual  rocket  exhausts.  The  reservoir  pressure  to  back 
pressure  relation  was  large,  as  high  as  350,  and  was  not  constant  but 
varied  markedly  with  time  during  an  experiment.  The  resulting  nozzle 
exit  pressure/back  pressure  ratio  varied  from  9.2  to  3.4.  iTie  tube 
configurations  considered  consisted  of  relatively  short  cylindrical 
tubes,  less  than  6.2  tube  internal  diameters  long,  with  a tube  radius 
to  nozzle  exit  radius,  rj./rg  of  1.2.  The  analytical  phase  of  the 
study  sought  to  verify  the  empirical  flow  model.  Based  on  the  data 
obtained  in  the  present  study  the  following  conclusions  were  made 
about  the  flow  field  in  the  cylindrical  tube: 

(1)  The  nozzle  exhaust  flow  along  the  entire  length  of  the 
tube  is  supersonic  with  the  wall  static  pressures  being 
greater  than  the  atmospheric  back  pressure.  The  base 
pressure  is  less  than  the  nozzle  exit  pressure  indi- 
cating that  underexpanded  nozzle  flow  exists. 

(2)  The  base  pressure/reservoir  pressure  relation  is  not 
constant  but  varies  for  fully  developed  supersonic  flow 
as  the  reservoir  pressure  varies.  The  ratio  of  base 
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pressure/reservoir  pressure  approaches  the  predicted 
static  equilibrium  value.  The  base  pressure  depends  on 
the  reservoir  pressure  and  the  resultant  flow  field 
characteristics  in  the  tube  are  dependent  upon  the  base 
pressure/reservoir  pressure  ratio. 

(3)  The  wall  pressure  measurements  indicate  an  abrupt 
pressure  rise  downstream  of  the  first  impingement  shock. 
For  the  conrltions  considered,  the  viscous  effects  have 
little  apparent  influence  on  the  impingement  shock 
structure  other  than  the  effects  on  base  pressure. 

(4)  The  flow  in  the  tube  is  Independent  of  tube  length. 

(5)  After  the  base  pressure  is  determined  from  an  existing 
theory,  the  location  of  the  first  impingement  shock  and 
the  wall  pressure  distribution  in  the  vicinity  of  the 
shock  wave  can  be  calculated  analytically, 

(6)  The  wall  pressure  peak  occurs  immediately  downstream 
of  the  first  impingement  shock  and  is  greatest  when 

the  reservoir  pressure  is  the  greatest.  Thus,  the 

1 

critical  stress  for  the  tube  design  can  be  readily 
calculated.  | 
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APPENDIX  A 


FINITE  DIFFERENCE  EQUATIONS  FOR  STEADY,  AXISYMMETRIC , 
POTENTIAL  FLOW  AND  STEADY,  AXISYMMETRIC, 
ROTATION  FLOW 


Equations  of  Characteristics 


The  following  equation  Is  valid  for  steady,  axlsymmetrlc, 
potential  flow: 


cot  4 


sin  u sin  0 dz 
cos(0  + |i)  r 


(A-1) 


For  steady,  axlsymmetrlc,  rotational  flow,  the  following  equation  Is 
valid: 


(A-2) 


Along  physical  characteristics  lines 


» tan  (e  ± n) 


(A-3) 


and 


dS  - 0 


Along  the  streamline  defined  by 


dz 


tan  e 


(A-4) 
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The  previously  mentioned  equations  are  derived  In  and  are 
taken  from  Ferrl  [26]. 

By  substituting  Equation  (A-3)  and  the  Prandtl-Meyar  function 
du  ■ cot  u ^ Into  Equations  (A-1)  and  (A-2) : 


d(„  +e)  . ^ 

sin(0  ± |i)  r 


(A-5) 


,,  -r  „ \ sin  u sin  6 dr  . . ao  ^ 

- ■si4"±"vt  t + 


(A-6) 


Since 


sin(0  ± 4)  =>  sin  6 cos  4 ± cos  6 sin  4 , 


sin  u cos  0 
sin(0  + 4) 


becomes 


cot  4 ± cot  0 


Equation  (A-5)  becomes 


d(v  + 0)  = 


cot  4 ± cot  0 r 


(A-7) 


and  Equation  (A-6)  becomes 


, . V 1 dr  , dS 

cot  n t cot  e T - ra 


(A-8) 


Steady,  Axlsymmetrlc,  Potential  Flow 
Finite  Difference  Equations 


Points  1 and  2 In  Figure  A-1  are  two  known  points  at  which 


the  flow  characteristics  are  known.  The  position  and  flow 


1 


(v  - 0),  - (v  - 0) 


2 1 


In 


^ ^cot  |i2  + ^2  *^3  '*'  ®s] 

(A-12) 


and 


(v  + 0),  - (v  + 0) 


11 


In-^  . 


|cot  lij^  - cot  + cot  (1^  “ cot  0^j  *4 

(A- 13) 


Equations  (A-10)  through  (A-13)  are  the  governing  character- 
istics equations  for  steady,  axLsymmetric,  potential  flow  in  finite 
difference  forji.  These  equations  are  identical  to  those  of  Owezarek 
[24]. 


Steady,  Axis )nnme trie,  Rotational  Flow 
Finite  Difference  Equations 


The 

pressure  drop 


relation  for  change  in  entropy  with  respect  to  total 
is 


<R 


(A- 14) 


Substitute  Equation  (A-14)  into  Equation  (A-6)  and  realize 
that  sin  2 p >•  2 sin  p cos  p gives 


d(v  + 0)  - 


cot  p,  ± cot  0 r 


dr  sin  2 u ^^0 

" 2 7 


(A-15) 
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Note  Equation  (A-3)  Is 


•jI  ■ tan(0  ± ji)  . (A-16) 

Because  the  rotational  flow  equations  will  be  used  behind  the 
shock  and  the  values  of  sotue  6's  and  r's  are  likely  to  be  zero,  the 
following  finite  difference  form  will  be  used  to  facilitate  numerical 
evaluation: 

AP  ^ 

A(v  + 0]- ^ (A-17) 

cot  |i  ± cot  0 r ' P 

g»tan[0±]I]  . (A- 18) 

The  finite  difference  form  of  Equations  (A-17)  and  (A-18)  are  the  same 
type  as  derived  by  Ruo  [27]. 

Equations  (A-17)  and  (A-18)  are  the  governing  characteristics 
equations  for  steady,  axlsymmetrlc,  rotational  flow  in  finite  differ- 
ence form. 

1 

■ ■/ 
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APPENDIX  B 


DETAILED  COMPUTATION  METHOD  FOR  AXIALLY  S'  MMETRIC 
ROTATIONAL  FLOW  REGION  BY  AN  APPROXIMATE  FINITE 
DIFFERENCE  METHOD  OF  CHARACTERISTICS 


Flow  Conditions  Behind  the  Impingement  Shock 

The  following  examples  of  the  use  of  the  rotational  flow 
equations  are  orientated  to  the  characterization  of  the  flow  behind 
the  Impingement  shock  wave.  Ir.  particular,  these  equations  in  con- 
junction with  the  equations  for  an  oblique  shock  will  be  used  to  deter- 
mine the  position  of  the  shock,  the  total  pressure  drop  across  the 
shock,  the  Mach  number,  and  flew  direction  behind  the  shock  [27,  28, 
29].  Figure  B-1  shows  the  angles  through  an  oblique  shock. 


WALL 


Figure  B-1.  Flow  Through  an  01)lique  Shock  Wave 


The  corditlons  behind  the  shock  at  the  jet  boundary-wall 
intersection  (position  I in  Figure  B-2)  and  the  shock  angle,  6,  are 
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readily  determined  from  the  oblique  shoc>  equations,  because  62,  the 

4 

flow  angle  behind  the  shock,  is  zero;  i.e.,  the  flow  turns  parallel 
with  the  tube  wall. 

Because  the  wall  constitutes  a streamline  for  Invlscld  flow, 
the  angle  62  ■ 0 and  P02/P0I  the  wall  Is  constant  (dS  « 0 along 

a streamilneK  However,  the  Mach  number  and,  consequently,  the  pres- 
sure along  the  wall  may  vary. 

Isentroplc  Flow  Characteristic  - Shock  !\*a ve  Intersect.i.on 

The  procedure  required  to  calculate  the  flow  down  stream  of 
the  shock  wave  at  point  III  Is  typical  of  the  procedure  for  points 
along  the  lmplugeme..t  shock  but  away  from  the  wall.  The  position  of 
point  III  can  be  obtained  from  geometry  because  the  shock  angle,  6,  Is 
known  at  point  II.  Flow  conditions  upstream  of  the  shock  at  point  III 
are  obtained  from  the  Isentroplc  flew  solution.  At  point  III,  a 
problem  arises  because  neither  the  Mach  number  nor  the  flow  angle 
behind  the  shock  Is  known.  The  oblique  shock  equations  require  that 
one  or  the  other  be  known.  A characteristics  net  must  be  constructed 
behind  the  shock  and  the  rotational  flow  and  oblique  shock  equations 
must  be  solved  simultaneously.  In  this  case,  they  are  solved  numeri- 
cally by  an  Iteration  process. 

Consider  Figure  B-2  where  point  III  is  close  to  the  wall  and 
the  flow  conditions  behind  the  shock  are  unknown.  For  illustration, 
assume  that  the  flow  coTiditions  along  the  waJl  ere  known  a function 
of  z out  to  B.  Hence,  if  the  position  of  C can  be  determined, 
the  flow  conditions  are  known.  Also,  assume  that  the  conditions  behind 
the  shock  at  point  II  are  known. 
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However,  the  conditions  behind  the  shock  at  point  IIX  are 
unknown.  Assume  that  62  at  point  III  is  equal  to  62  at  point  II  and 
solve  the  shock  equations  to  determine  the  flow  conditions  behind  the 
shock. 

The  flow  conditions  must  also  satisfy  the  rotational  flow 
characteristics  net  behind  the  shock.  Solve  the  characteristics 
equations  along  the  characteristic  C-III  in  terms  of  by  the  finite 

difference  equations: 


[v  + eiiji  - ( vl^ 


cot 


\\ 

sin  [ix^  + ^0 


Cot 


III 


'■>111  \ 


‘0. 


P P 

P P 

^0  0 


and 


*^111  ■ ’^t  ^ pin  f^IIll 

__ . j^_  _j 


Pin  ~ ^t] 

•|pni+ 


(B-l) 


(B-2) 


These  equations  cannot  be  solved  because  the  characteristic- 
wall  intersection  position,  , is  unknown.  If  z^  were  known,  the 
flow  conditions  at  C would  be  known  (based  on  an  original  assumption). 
For  the  first  approximation,  let  the  flow  conditions  at  C equal  to  the 
flow  conditions  of  the  previous  characteristic-wall  intersection  A, 
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i.e.,  and  By  solving  Equation  (B-2)  , a new  z-position 

for  point  C will  be  generated.  Because  the  flow  conditions  along  the 
wail  are  assuned  a function  of  z,  the  flow  conditions  are  also  known. 
The  new  values  are  substituted  Into  Equation  (B-2).  This  process  Is 
repeated  until  ^ c is  a small  acceptable 


value . 


Note  the  first  approximation  of  is  obtained  using  the 

Mach  angle,  found  by  the  oblique  shock  equations.  The  next 

approximation  for  Vjjj  is  found  by  using  the  Mach  angle  corresponding 
to  the  previous  Prandtl-Meyer  function.  Note  that  the  position  a'.  C 
also  changes,  so  the  iteration  procedure  for  position  C must  be 
conducted. 


This  procedure  is  repeated  until  - (v,,,)  < c. 

m (»-l) 

Equation  (B-1)  is  solved. 

Now,  if  both  conditions  are  satisfied  (oblique  shock 
equations-rotational  flow  characteristics  net),  Vg  g “ ''m  C ‘ 
not,  02  is  Incremented  by  a small  value,  A62,  with  the  attendant  iter- 
ations until  this  condition  is  met. 

Rotational  Characteristic  - Tube  Wall  Intersection 


The  conditions  at  D on  the  wall  must  be  satisfied  so  that  the 
solution  may  progress  along  the  wall  for  points  past  B.  This  is  done 
by  solving  the  finite  difference  equations^  for  the  characteristics 
extending  from  point  III  to  D. 


r 


[vId  - (V  - eljii 


cot  2 


‘III 


+ cot 


III 


- "ill] 

2 [^'t  *‘lll] 


ain  [mi^I  + Up]  [ ^0  **o  J 


'III 


P P 

1 1 °i  , °m 
'o  ‘■o 


(B-3) 


^t  ~ *^III 
*D  ” ^III 


tan 


pill  ^^III 
[ 2 2 


(B-4) 


To  solve  these  equations,  Vq  must  be  known.  In  this  case, 
assume  for  the  first  estimate  Vp  • Vg  and  Iterate  until 


[''D(m-l)  ■ 'D(m)j  ^ ^ 

An  equation  must  be  used  to  describe  the  variation  of  Mach 
nun^er  along  the  wall  as  a function  of  z. 

Suppose  the  z position  of  C Is  known,  but  v^.  Is  unknown. 
Remember  the  flow  conditions  and  positions  of  A and  B are  known.  Now, 
If  A and  B are  sufficiently  close,  v^,  may  be  obtained  by 


'’C  “ "a  + 


V •»  V 

. B / X 

'c  ■''A'^vr^('c.-^A) 


(B-5) 
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I n Itlal  Rotational  Characteristic  - Shock  Wave  Intersection 


In  Che  solution  for  polnc  III  behind  the  shock  and  D on  Che 
wall,  Che  flow  conditions  along  Che  wall  were  known  out  to  B along  Che 
wall.  However,  in  solving  for  point  II,  only  Che  conditions  behind  Che 
shock  at  point  1 and  the  vail  streamline  are  known. 

This  problem  is  solved  numerically  by  letting  the  conditions 
behind  Che  shock  at  I equal  to  the  conditions  at  A (first  approxi- 
mation) and  solving  for  Che  conditions  behind  Che  shock.  We  then 
project  to  Che  wall  at  B through  Equations  (B-4)  and  (B-5) . Now 
Equation  (3-5)  is  solved  to  obtain  a better  estimate  of  v^.  This 
process  is  repeated  tmcll 

pACm-l)  ■ V(in)]  ^ ^ 

Accuracy  Discussion 

The  primary  purpose  of  this  exercise  is  to  determine  the 
pressure  along  Che  wall  aft  of  Che  impingement  shock.  Prior  to  this 
time,  approximate  methods,  such  as  Bauer  [17]  simplified  theory, 
assumed  that  Che  flow  behind  Che  shock  was  two-dimensional,  was 
deflected  parallel  Co  the  wall,  and  that  Che  pressures  behind  Che 
shock  were  lirpressed  on  Che  wall.  He  bases  the  Justification  for  this 
simplified  theory  on  the  premise  that  the  variation  of  the  flow  angle 
with  respect  to  that  assumed  is  small,  l.e.,  ^2  * 

neglected  second-order  terms  6*’.  Bauer  also  limits  the  applicability 
of  this  theory  to  the  region  near  the  shock. 
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A complete  computer  program  which  would  define  the  flow  via 
the  method  of  characteriatics  for  multiple  shock  flows  was  Initiated 
by  Ruo  [27].  However,  Ruo  was  unable  to  complete  this  formidable  task 
because  of  a lack  of  funding.  A more  exact  techr  '.que  is  well  within 
the  state-of-the-art  of  computer  science,  provided  one  has  the  time 
and  funds. 

The  parameter  in  the  equation  which  differentiates  axisym- 
metric  flow  from  two-dimensional  plane  flow  is 

1 ^ 

cot  4 ± cot  0 r 

Note  as  the  flow  angle,  8,  approaches  zero,  the  cot  6 
approaches  infinity,  and  if  the  other  terms  remain  finite,  the  overall 
term  tends  toward  zero.  Now,  if  the  point  under  consideration  is  near 
the  wall,  dr  or  Ar  is  small  and  r is  large,  and  6 1?  small,  this  term 
can  be  neglected  and  the  flow  reverts  to  two-dimensional  plane  flow  as 
was  Bauer's  assumption.  However,  in  cur  case,  6 is  small  but  maybe  not 
that  small,  and  dr  or  Ar  can  be  large,  depending  on  the  distance  of 
the  point  under  consideration  from  the  wall  streamline.  Tlie  question 
is  how  large  can  the  relation  be,  or  how  close  to  the  tube  axis 

will  the  calculation  yield  acceptable  answers . Certainly  the  accept- 
ability depends  on  0 which  is  not  known.  Another  problem  exists,  con- 
sidering that  APq  may  not  be  small  even  if  dr  is.  If  APq  is  large, 
then  acceptable  results  are  obtained  only  if  the  entropy  term  can  be 
Integrated,  l.e.. 
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sin  2 u 

J 2, 


To  do  this,  sin  2 u must  be  a constant  or  not  a function  of  Fq.  This, 
In  a way.  Is  a contradiction,  because  In  the  axlsymmetrlc  term,  6 
appears  with  p.  So,  If  p tends 'toward  a constant,  then  this  tern  can- 
not be  Integrated,  which  Is  required  If  Ar  Is  large.  However,  If  In 
reality  p Is  small  and  p does  not  vary  too  strongly  along  a charac- 
teristic, this  term  will  tend  toward  a constant  and  reasonable  esti- 
mates for  large  values  of  At  and  APq  may  be  obtained.  Because  a shod’, 
tends  to  process  a nonuniform  flow  field  toward  these  requirements, 
the  prospect  seems  good  that  reasonable  estimates  will  be  obtained. 
However,  in  reality,  a comparison  of  the  analytical  estimates  to 
experimental  values  seems  to  be  the  best  method  of  determining  an 
accuracy  limitation  which  would  answer  the  basic  question. 
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APPENIdX  C 


FLOW  CHAKACTE VIST ICE 


* t. 

Tat.le  C-I.  Cold  Flow  Isentroplc  Characteristics 0.02A4,  Y ■ 1.4 

pq 

(For  Characteristics  Net  in  Figure  C-1) 


Points 

Mach 

No. 

Flow 

Angle 

(dag) 

P 

PO 

Poin^" 

Kach 

No. 

H 

■ 

P 

Po 

00 

2.777 

20.000 

0.03815 

22 

3.031 

21.971 

0.02598 

01 

2.848 

21.470 

0.03420 

23 

3.114 

23.454 

0.02297 

02 

2.921 

22.940 

0.02066 

24 

3.199 

24.936 

0.02024 

03 

2.996 

24.410 

0. 02739 

25 

3.132 

23.156 

0.02237 

04* 

3.073 

25.880 

0.02440 

i 

26* 

3.073 

21.527 

0.02440 

10 

2.825 

19.500 

0.03545 

30 

2.924 

18.500 

0.03052 

11 

2.899 

20.978 

0.03168 

31 

3.005 

19.997 

0.02703 

12 

2.976 

22.455 

0.02824 

32 

3.088 

21.488 

0.02386 

13 

3.054 

23.931 

0.02510 

33 

3.175 

22.978 

0.02099 

14 

3.136 

25.407 

0.02224 

34 

3.265 

24.465 

0.01839 

15* 

3.073 

23.641 

0.02|440 

35 

3.192 

22.672 

0.02047 

1 

36 

3.129 

21.031 

0.02247 

20 

2.874 

19.000 

0.03291 

37* 

3.073 

19.519 

0.02440 

21 

2.951 

20.487 

0.02928 

* Jet  Boundary 


jHKrntmm  ■ umi 


PRECEDIN3  PASS 

— < .. , ■— 
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Table  C-I.  Continued 


Flow 

Flow 

Mach 

Angle 

P 

Mach 

Angle 

P 

Points 

No. 

(deg) 

Po 

Points 

No. 

(deg) 

Po 

40 

2.975 

18.000 

0.02827 

60 

3.079 

17.000 

0.02418 

41 

3.059 

19.507 

0.02492 

61 

3.172 

18.530 

0.02109 

42 

3.147 

21.006 

0.02189 

62 

3.267 

20.043 

0.01833 

43 

3.114 

22.502 

0.01914 

63 

3.368 

21.554 

0.01584 

44 

3.332 

23.996 

0.01667 

64 

3.473 

23.060 

0.01363 

45 

3.253 

22.189 

0.01871 

65 

3.381 

21.225 

0.01554 

46 

3.185 

20.335 

0.02067 

66 

3.303 

19.544 

0.01741 

47 

3.126 

19.013 

0.02.'>56 

67** 

3.254 

18.437 

0.01868 

Jet 

Bound- 

ary-Wall 

Intersec- 

70 

3.133 

16.500 

0.02231 

tion** 

3.073 

18.197 

0.02440 

71 

3.230 

18.049 

0.01935 

48 

3.096 

18.216 

0.02359 

72 

3.330 

19.563 

0.01673 

73 

3.435 

21.082 

0.01438 

50 

3.027 

17.500 

0.02616 

74 

3.546 

22.594 

0.01229 

51 

3.115 

19.018 

0.02294 

75 

3.447 

20.744 

0.01414 

52 

3.206 

20.524 

0.02004 

76 

3.363 

19.049 

0.01594 

53 

3.302 

22.028 

0.01743 

77** 

3.351 

18.792 

0.01622 

54 

3.402 

23.528 

0.01509 

55 

3.316 

21.707 

0.01707 

80 

3.189 

16.000 

0.02057 

56 

3.243 

20.039 

0.01898 

81 

3.289 

17.552 

0.01775 

57 

3.180 

18.507 

0.02083 

82 

3.394 

19.084 

0.01525 

58** 

3.171 

18.268 

0.02112 

83 

3.505 

20.610 

6.01302 

**  Shock  Intersection 


/■ 


/ 

/' 
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Table  C-I 


Points 

Mach 

Flow 

Angle 

(deg) 

P 

No. 

?0 

84 

3.621 

22.128 

0.01105 

85 

3.515 

20.283 

0.01283 

86** 

3.480 

19.600 

0.01349 

90 

3.245 

15.500 

0.01894 

91 

3.350 

1/.062 

0.01625 

92 

3.460 

18.606 

0.01387 

93 

3.576 

20.140 

0.01177 

94** 

3.696 

21.632 

0.00995 

100 

3.303 

15.000 

0.01741 

101 

3.413 

16.573 

0.01485 

102 

3.528 

18.128 

0.01260 

103 

3.650 

19.670 

0.010G2 

104** 

3.656 

19.743 

0.01053 

**  Shock  Intersection 


S 


Points 

Mach 

No. 

Flow 

Angle 

(deg) 

P 

^0 

110 

3.361 

14.500 

0.01599 

111 

3.476 

16.084 

0.01356 

112 

3.598 

17.651 

0.01142 

113** 

3.628 

18.022 

0.01095 

120 

3.422 

14.000 

0.01466 

121 

3.542 

15.597 

0.01235 

122** 

3.610 

16.443 

0.01123 

130 

3.483 

13.500 

0.01343 

131** 

3.600 

14.995 

0.01138 

140 

3.5.46 

13.000 

0.01228 

141** 

3.595 

13.607 

0.01146 

150** 


3.604 


12.500 


0.01132 


{**/*)  NOIllSOd  1VIOVH 
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JET  BOUNDARY 


TUBE  WALL 


AXIAL  POSITION  FROM  NOZZLE  EXIT  (|z  - zjl  /r^] 

Figure  C-3 . Concluded 
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p 

Table  C-II.  Hot  Flow  Isentropic  Characteristics —=  0.0118,  Y • 1.2 

0 

(For  Characteristics  Net  in  Figure  C~1Y 


Points 

Mach 

No. 

Flow 

Angle 

(deg) 

P 

^*0 

Points 

No. 

- 

Flow 

Angle 

(deg) 

P 

"o 

00 

2.941 

20.000 

0.02493 

30 

3.053 

18.500 

0.0203 

01 

3.039 

22.627 

0.02084 

31 

3.164 

21.183 

0.0166. 

02 

3.141 

25.254 

0.01734 

32 

3.278 

23.859 

0.0135. 

03 

3.246 

27.880 

0.01434 

33 

3.398 

26.530 

0.0109( 

04* 

3.354 

30.507 

0.01180 

34 

3.524 

29.194 

0.0087( 

35 

3.457 

27.064 

0.0098f 

10 

2.978 

19.500 

0.02330 

36 

3.402 

25.168 

0.0108; 

11 

3.080 

22.144 

0.01935 

37** 

3.354 

24.219 

0.0118( 

12 

3.186 

24.787 

0.01598 

13 

3.296 

27.429 

0.01311 

40 

3.092 

18.000 

0.01895 

14 

3.410 

30.068 

0.01068 

41 

3.206 

20.703 

0.01541 

15* 

3.354 

27.981 

0.01180 

42 

3.326 

23.396 

0.01242 

43 

3.451 

26.081 

0.00992 

20 

3.015 

19.000 

0.02177 

44 

3.583 

28.757 

0.00784 

21 

3.122 

21.664 

0.01796 

\ 

45 

3.510 

26.607 

0.00892 

22 

3.232, 

24.323 

0.01471 

1 

46** 

3.451 

24.711 

0.00992 

23 

3.347 

26.979 

0.01196 

. 

24 

3.466  ' 

29.631 

0.00964 

50 

3.131 

17.500 

0.01766 

25 

3.405 

1 27.523 

0.01129 

51 

3.249 

20.224 

0.01425 

26* 

3.354 

(25-644 

.-i, 

0.01180 

52 

3.374 

22.935 

0,01139 

* Jet  Boundary 
**  Shock  Intersection 
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Table  C-II.  Concluded 


IKSi 

Mach 

P 

Mach 

Points 

No. 

PO 

Points 

No. 

53 

3.504 

24.634 

0.00901 

90 

3,291 

54 

3.643 

28.320 

0.00705 

91 

3.428 

55 

3.564 

26.150 

0.00810 

92 

3.573 

56 ‘V 

3.540 

25.437 

0.00846 

93** 

3.594 

60 

3.170 

17.000 

0.01644 

100 

3.333 

61 

3.293 

19.745 

0.01317 

101 

3.475 

62 

3.422 

22.474 

0.01044 

102** 

3.569 

63 

3.559 

24.187 

0.00818 

110 

3.375 

64 

3.704 

27.884 

0.00633 

111 

3.522 

65** 

3.663 

26.069 

0.00679 

112** 

3.554 

70 

3.210 

16.500 

0.01530 

71 

3.337 

19.267 

0.01216 

120 

3.417 

72 

3.472 

22.014 

0.00955 

121** 

3.545 

73 

3.614 

24.7^2 

0.00741 

130 

3.460 

74** 

3.677 

25.873 

0.00663 

131** 

3.539 

oO 

3.250 

16.000 

0.01422 

140 

3.504 

81 

3.38.' 

18.790 

0.01121 

141** 

3.540 

82 

3.522 

21.555 

0.00873 

83** 

3.628 

23.521 

0.00723 

150** 

3.546 

15.500 

0.01321 

18,314 

0.01033 

21.097 

0.00797 

21.474 

i 

0.00768 

15.000 

0.01226 

17.839 

0.00951 

19.611 

0.00803 

14.500 

0.01137 

17.364 

0.00874 

17.956 

0.00825 

14.000 

0.01054 

16.434 

0.00838 

13.500 

0.00975 

14.977 

0.00847 

13.000 

0.00902 

1^.055 

0.00847 

12.500 

0.00836 

* Jet  Boundary 
**  Shock  Intersection 


107 


108 


gure  C-2.  IsenCroplc  Characteristics 


Figure  C-2.  Concluded 
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